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Abstract 

We study theoretically the role of ageing in the rheology of soft materials. 
We define several generalized rheological response functions suited to ageing 
samples (in which time translation invariance is lost). These are then used to 
study ageing effects within a simple scalar model (the "soft glassy rheology" 
or SGR model) whose constitutive equations relate shear stress to shear 
strain among a set of elastic elements, with distributed yield thresholds, 
undergoing activated dynamics governed by a "noise temperature" , x. (Be- 
tween yields, each element follows affinely the applied shear.) For f < x < 2 
there is a power-law fluid regime in which transients occur, but no ageing. 
For x < 1, the model has a macroscopic yield stress. So long as this yield 
stress is not exceeded, ageing occurs, with a sample's apparent relaxation 
time being of order its own age. The (age-dependent) linear viscoelastic loss 
modulus G"(u>,t) rises as frequency is lowered, but falls with age t, so as to 
always remain less than G'(to,t) (which is nearly constant). Significant age- 
ing is also predicted for the stress overshoot in nonlinear shear startup and 
for the creep compliance. Though obviously oversimplified, the SGR model 
may provide a valuable paradigm for the experimental and theoretical study 
of rheological ageing phenomena in soft solids. 



1 Introduction 



Many soft materials, such as foams, dense emulsions, pastes and slurries, display 
intriguing features in their low frequency shear rheology. In oscillatory shear, for 
example, their viscoelastic storage and loss moduli, G'{uj) and G"(lu), are often 



weak power laws of shear frequency (Mackley et al., 1994|; Ketz et al., 1988 


; Khan 


et al, 1988; Mason et al., 1995|; Panizza et al., 1996 


; Hoffmann and Rauscher, 
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1993; Mason and Weitz, 1995), while their nonlinear stress response a to shear 
strain of constant rate 7 is often fit to the form a = A + Bj n (known as the 
Herschel-Bulkley equation, or when A = 0, the power-law fluid) QHoldsworth 



1993; Dickinson, 1992; Barnes et al,, 1989| ) . The fact that such a broad family 
of soft materials exhibits similar rheological anomalies is suggestive of a common 
cause, and it has been argued that these anomalies are symptomatic of the generic 
presence in such materials of slow, glassy dynamics ( jSollich et al., 1997; Sollich 



1998| ), Indeed, all the above materials share features of structural disorder and 



metastability: large energy barriers impede reorganization into states of lower free 
energy because this would require rearrangement of local structural units, such 
as the droplets in a dense emulsion. The term "soft glassy materials" (SGM's) 
has been proposed to describe such materials ( [Sollich et al., 19"97; Sollich, 1998| ). 

Glassy dynamics are often studied using hopping (trap) models, in which sin- 
gle particle degrees of freedom hop by an activated dynamics, in an uncorrelated 
manner, through a random free energy landscape ( Bouchaud, 1992; Monthus and 
Bouchaud, 1996[ ). By incorporating strain degrees of freedom into such a de- 
scription, Sollich and coworkers ( pollich et al., 1997|; Sollich, 1998| ) proposed a 
minimal model, called the "soft glassy rheology" (SGR) model, which appears to 
capture several of the rheological properties of SGM's, although (for simplicity) all 
the tensorial aspects of viscoelasticity are discarded. The model exhibits various 
regimes depending on a parameter x (discussed in more detail below) representing 
the "effective temperature" for the hopping process. When this is small (x < 1) 
the model exhibits a glass phase which shows some interesting properties above 
and beyond the power-law anomalies in viscoelasticity mentioned above. Specif- 
ically, the model shows ageing behaviour, its properties depend on the elapsed 
time since a sample was prepared. This is because the population of traps visited 
never achieves a steady state; as time goes by, deeper and deeper traps dominate 
the behaviour (a phenomenon known as "weak ergodicity breaking"). Broadly 
speaking, the system behaves as though its longest relaxation time is of order its 
own age. 

The success of the SGR model in accounting for some of the generic flow prop- 
erties of SGM's suggests that a detailed investigation of its ageing behaviour, and 
the effect this has on rheology, is now worthwhile. Ageing has been intensively 
studied in the context of spin glasses (|Bouchaud and Dean, 1995; Cugliandolo| 
and Kurchan, 1995; Bouchaud et al., 199 j ), although some of the earliest experi- 
mental investigations of it involved rheological studies of glassy polymers ( [Struik 



1978| ). But we know of no previous theoretical work that explores the link be- 



tween ageing phenomena and rheological properties within an explicit constitutive 
model. A particular added motivation is that detailed experiments on rheolog- 
ical ageing, in a dense microgel suspension, are now underway flCloitre, 1999 ). 
Although various kinds of ageing effects are often observable experimentally in 
soft materials, they have rarely been reported in detail. Instead they tend to be 
regarded as unwanted obstacles to observing the "real" behaviour of the system, 
and not in themselves worthy of study. But this may be illusory: ageing, when 
present, can form an integral part of a sample's rheological response. For exam- 
ple, the literature contains many reports of viscoelastic spectra in which the loss 
modulus G"{uj), while remaining less than the (almost constant) storage modulus 
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G'{uj) in a measured frequency window, appears to be increasing as frequency is 
lowered (see Fig. HJ). The usual explanation ( Kossuth et al., 1999 ) is that some 
unspecified relaxation process is occurring at a lower frequency still, giving a loss 
peak (dashed), whose true nature could be elucidated if only the frequency win- 
dow was extended. This may often be the case, but an alternative explanation, 
based on our explicit calculations for the SGR model, is shown by the thin solid 
lines. No oscillatory measurement can probe a frequency far below the reciprocal 
of the sample's age; yet in ageing materials, it is the age itself which sets the 
relaxation time of whatever slow relaxations are present. Accordingly, the puta- 
tive loss "peak" can never be observed and is, in fact, a complete figment of the 
imagination. Instead, a rising curve in G"{uo) at low frequencies will always be 
seen, but with an amplitude that decreases as the system gets older (typically 
ensuring that G"(u) never exceeds G'(oS)). Of course, we do not argue that all 
published spectra resembling those of Fig. || should be interpreted in this way; 
but we believe that many should be. The widespread reluctance to acknowledge 
the role of ageing effects in much of the rheological literature suggests that a full 
discussion of these could now be valuable. An exception has been in the literature 
on ageing in polymeric glasses, especially the monog raph by |Struik (19781) : we 
return shortly to a brief comparison between that work and ours. 

The SGR model is simple enough to allow a fairly full exploration of the link 
between ageing and rheology. As well as providing some quantitative predictions 
of rheological ageing, this allows a broader discussion of the conceptual frame- 
work within which rheological data for ageing systems should be analysed and 
interpreted. This conceptual framework is broader than the SGR model itself; 
for example it is known that ageing concepts developed for spin-glass dynamics 
can also be applied to problems of domain-growth and coarsening ( pouchaud 
et al., 1998|) . Many soft solids, such as polydomain defect textures in ordered 
mesophases of copolymers or surfactants, may show ageing through such coars- 
ening dynamics, or through glassy rearrangement of domains, or both. While the 
SGR model is intended to address only the second feature, the broader conceptual 
framework we present can allow for both mechanisms (in which case a superposi- 
tion of ageing dynamics with different timescales may result; see Eq. (|2|) below). 

Thus we begin in Sees. |2| and ||| by briefly introducing rheology and ageing 
respectively. Then in Sec. || we review the SGR model, and discuss the origin of 
its glass transition and the nature of the glass phase. We also briefly describe its 
rheology under non-ageing conditions; this is discussed fully elsewhere ( [Sollich 
et al., 1997 ; Sollich, 1998| ). In Sec. || we give a general discussion of ageing 
within the SGR model, which sets the stage for our new results for the linear 
and nonlinear rheological response of the SGR model in regimes where ageing 
cannot be neglected. The results for controlled strain conditions are presented 
and discussed in Sec. ||; those for controlled stress, in Sec. [/]. We close in Sec. || 
with a brief summary and our conclusions. 

We now discuss the connection between our work and that of Struik (1978| ) on 
polymeric glasses. Struik presented many experimental results for such systems, 
and gave a coherent qualitative explanation of their ageing in terms of a slow 
relaxation of the free volume in the system below the glass point. However, 
he did not propose a model constitutive equation for this or any other class of 
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Figure 1: Sketch of ageing scenario for dynamic moduli G' (top) and G" (bot- 
tom). Assuming that data shown by the thick solid lines has been measured, the 
conventional interpretation would be that G" has a loss peak at low frequencies 
oj (dashed) outside the experimentally accessible range. But the rise of G" to- 
wards small uj could also be due to ageing: If the experiment is repeated at later 
times, the thin solid lines (calculated from the SGR model at noise temperature 
x = 0.7; see Sec. 3.1.2 ) might be measured; the putative loss peak would then 
always be unobservable because no measurement is possible at frequencies below 
the inverse age of the system. 



ageing material. He argued that the effective relaxation time of a system of age 
t w (the 'waiting time' since sample preparation) varies r(t w ) = i^T 1_M , where 
To is a microscopic time and \i ~ 1; apart from this uniform rescaling of (all) 
the rheological relaxation time(s), the material properties are almost invariant 
in time. (This is his 'time waiting-time superposition' principle; we show below 
that the SGR model offers a concrete example of it, with fj, = 1.) We do not 
expect the SGR model, which makes no mention of the free-volume concept, 
to be particularly relevant to polymeric glasses; nonetheless, various points of 
contact with Struik's work are indicated below. 



2 Rheology 

Here we review the basic principles of rheology. Unlike most in the literature, 
our formulation does not assume time translational invariance (TTI); parts of 
it may therefore be unfamiliar, even to rheologists. The formalism allows in 
principle an arbitrary dependence of the material properties on time; we defer 



to Sec. H a discussion of what more specific form this dependence might take in 
materials which exhibit actual ageing effects (rather than other, more trivial time 
dependencies). 



2.1 Constitutive Properties 

Rheology is the study of the deformation and flow properties of materials. In 
general, deformation can comprise volume changes, extensional strain, and shear 
strain; here we consider incompressible materials and assume that only shear 
strains arise. A system's shear stress a(t) then depends functionally on its strain 
rate history j(t' < t), where 7 is the strain rate. Conversely, 7(i) can be expressed 
as a functional of the preceding stress history. A specification of either type is 
referred to as a constitutive equation. In general, of course, the constitutive 



equation is a relationship between stress and strain tensors; see Doi and Edwards 



(1986| ) for an introduction. We ignore the tensorial aspects here, because the 



model we describe later is too simple to address them. 



2.2 Step Strain 

A standard rheological test consists of suddenly straining a previously undeformed 
material by an amount 70. Suppose this to be done at time t w : then ^(t) = 
7o©(t — i w )j where G is the usual step function. (For the moment, i w is an 
arbitrary time label, but later we will take it as the time that the strain is 
applied, relative to the preparation of the sample in some prescribed state, at 
time zero.) In general the response can be written 

a(t) = l G(t - t w ,t w ; 70) (1) 

thereby defining the step strain response, G(t—t w ,t w ; 70). Note that, by causality, 
G vanishes for negative values of its first argument. 



2.3 Time Translation Invariance; Linearity 

If the material properties of a sample have TTI, then the time t w of the initial 
step strain is irrelevant; the response function G{t — t w ,t w ;~/Q) can be written 
G(t — t w ;7o) and depends only on the elapsed time since the step strain was 
imposed. 

It is particularly important to recognize that TTI is a quite separate issue 
from the linearity of the material's response to stress. Even when TTI is absent, 
in the small deformation limit (70 — > 0) , a regime may exist for which a is linearly 
related to 70: 

lim G(t-t w ,t w ; 70) = G(t-t w ,t w ) (2) 
70^0 

The system's stress response is then linearly proportional to strain amplitude (in 
the sense that doubling the strain at all earlier times will cause the stress to be 
doubled), even if it also depends on (say) the sample's age relative to an absolute 
time of preparation. 
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Only by assuming both linearity and TTI do we obtain 



a(t) = 7o G(t - t w ) (3) 

where the function G(t) is called the time-dependent modulus, or the linear stress 
relaxation function, of the material. If a linear material with TTI is subjected to 
a small time-dependent strain j(t), then by decomposing this into a sequence of 
infinitesimal step strains, one finds 

a(t)= C G(t-t')j(t')dt' (4) 



which is, for a linear material with TTI, the constitutive equation between stress 
and strain. In the steady state (i.e., for constant strain rate 7) one recovers: 

POO 

<r = 7 / G{t")dt" (5) 



The integral, whenever it exists, defines the material's zero-shear viscosity rj. For 
many soft materials, however, G(t) decays to zero so slowly that the integral 
diverges. In this case, there can be no regime of linear response in steady shear 
flow, although there may be a linear regime in, say, oscillatory shear. 

Note that there is no unique extension of (|j) to the nonlinear case; only in 
the linear regime can one superpose the contributions from each small strain 
increment in this fashion. (In some models, the stress for a general flow is indeed 
written as an integral involving the nonlinear step strain response ( [Bernstein 



et al., 1963); but this is not generally valid.) On the other hand, (|j) is easily 



extended to the case where TTI is absent: 

a{t) = [ G(t-t',t')j(t')dt' (6) 



which represents the most general form of a (nontensorial) linearized constitutive 
equation. 



2.4 Behaviour of the Linear Response Function 

The principle of causality demands that the response function G(t — t w ,t w ) is 
zero for times t < t w . At t = i w , when strain is applied, G typically increases 
very rapidly (in effect discontinuously) to a value Go which represents an instan- 
taneous elastic response with modulus Go. Thereafter, G(t — t w ,t w ) is (almost 
always) a decaying function of its first argument: the more nearly the material 
approximates a viscous liquid, the more rapidly will the stress decay. Specializing 
to the TTI case, we recall that for a purely Newtonian liquid of viscosity rj, the 
function G(t) approaches a delta function r]5(t). (This shows that Go can be 
infinite so long as the subsequent decay is rapid enough.) On the other hand an 
ideally Hookean elastic material has G(t) = Go, the static shear modulus: in this 
case the induced stress will never decay. (Note that properly one should write 
G(t) = Go©(t); the extra factor of O(t), implied by causality, is omitted here and 
below.) 
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Both the Newtonian fluid and the Hookean solid are idealized limiting cases; 
most real materials display behaviour intermediate between these limits and are, 
on some timescale at least, viscoelastic. For the soft materials of interest to us, the 
relevant timescale is readily observable in rheological experiments. The simplest 
(TTI) example of viscoelasticity is the Maxwell fluid, which is solid-like at short 
times and liquid at longer ones, with a simple exponential response function 
G(t) = Goexp(— t/r) connecting the two. (Its viscosity obeys r\ = Got.) This 
behaviour is seen in a few experimental systems ( Cates and Candau, 199C| ), but 



more often one has G(t) = Go/j,(t) where the memory function fj,(t) is not a single 
exponential. In many materials it is possible to identify a longest relaxation time 
via t~I x = — lim^oo log /z(i)/t. However, in several important cases, such as 
a pure power law relaxation, fj,(t) ~ t~ y , the required limit does not exist; the 
longest relaxation time is infinite. 

2.5 Creep Compliance 

Arguing along parallel lines to those developed above, one can in general write 
the strain response to a step stress a(t) = ao<d(t — i w ) as 

7(i) = a J(t - t w ,t w ;a ) (7) 

The linear creep compliance J(t — t w , t w ) is then found by letting gq — > (assum- 
ing this limit exists). This is the main rheological function considered by jBtruik 



1971 ) and the one most relevant to studies of ageing in polymeric glasses (since 



these are often used, for example, as structural components subject to steady 
loads). 

In the presence of TTI the linear compliance reduces to a function of one 
time variable, J[t — t w ). For the examples of a viscous liquid, an elastic solid, 
and a Maxwell material we have (again omitting factors of @(t)) J(t) = t/r), 
J(t) = 1/Gq, and J(t) = 1/Gq + t/r), respectively. For any material with TTI, 
the zero-shear viscosity r/ is defined experimentally as the limiting ratio of stress 
to strain rate long after application of an infinitesimal step stress; it therefore 
obeys r/ _1 = lim^oo dJ(t)/dt, which may be shownQ to be equivalent to (|5|). A 
finite viscosity, requires, of course, that the limit is finite; this is discussed further 
in Sec. below. 

2.6 Viscoelastic Spectra 

A common experiment is to apply a steady oscillatory strain and measure the 
resulting stress, or vice versa. For example, the choice 

7 (t) = 6(i - t s )Re [ 7 oe i( ^ } ] (8) 

describes an oscillatory flow started at time t s and continued up to (at least) the 
time t at which the stress is measured. Using the linear constitutive equation for 



1 The given limit may also be written r\ 1 = lim^^o icoJ* (lj) which, by reciprocity of J* and 



G* , implies r\ — lim^^o G* {u>)/iu). The last definition is equivalent to (^). See Sec. 2.6 
definitions of J* and G* . 



for 
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a system with TTI (||), we have 



a(t) = Re ' 



Re 



l0 iuj [ e i{M) G{t-t')dt' + lQ e i ^ +ujts) G{t-t i 

Jt s 



(9) 



where the second term accounts for any step strain arising at the switch on time 
t B . As the number of cycles becomes very large (uj(t — t s ) 3> 1), transient effects 
become negligible, and the stress settles to a simple oscillatory function of time. 
In this steady-state limit we can write a(t) = Re [G*(u))j(t)] whereQ 

POO 

G*{uj) =iuj / e- iu3i G(t)dt (10) 
Jo 

which is, to within a factor iui, the Fourier transform of the stress relaxation 
modulus G(t). Traditionally one writes 

G*(u) = G'(uj) + iG"(uj) (11) 

where G', G" are called respectively the storage and loss moduli of the material, 
and measure the in-phase (elastic) and out-of-phase (dissipative) response to an 
applied strainP]. 

Clearly one can reach an identical steady state by applying a small amplitude 
oscillatory stress and measuring the resulting strain. This defines a new response 
function J*(u>) via j(t) = Re [J* (to) a (t)], which is evidently just the reciprocal 
of G*(lo). But by an argument similar to that given above for ( JlOj ) one also 
has J*(uj) = iio Jq 00 e~ lu)t J(t) dt. Hence, within the linear response regime of a 
system with TTI, knowledge of any one of G(t), J(t),G* (lo), J*(u>) is enough to 
determine the other three. (Of course, this ignores any practical limitations on 
the time and frequency domains accessible by experiment.) 

Beyond the linear response regime, it is sometimes useful to define G*{uj\^q) 
and J*(lo; <7o) from the response to a finite amplitude oscillatory shear. However, 
the interest in these quantities is more limited since, whenever the strain depen- 
dence is nontrivial, there is no analogue of ( |i~0| ) relating the nonlinear oscillatory 
response to that in step strain or step stress. 



2.7 Viscoelastic Spectra without TTI 

The proper definition of linear viscoelastic spectra for systems without TTI is 
more subtle, and is to some extent a matter of choice. Let us envisage again the 
following idealized experiment: (i) the sample is prepared in a known state at 

2 If G(t) has a non-decaying contribution G(t — » oo) > 0, the form G*(lj) — G(0) + 
J" °° e~ lut G'(t) dt, derived from Q by integration by parts, should be used instead of (p"o|). 
The same relation can be obtained from ( ]lo| ) by inserting a regularizing factor e~ tt and taking 
the limit e — > 0. This corresponds to an oscillatory strain that is switched on by very slowly 
increasing its amplitude. 

3 Many commercial rheometers are configured to deliver the storage and loss spectra auto- 
matically, from a measurement of the amplitude and phase relations between stress and strain 
in steady state. 



s 



time zero; (ii) a small amplitude oscillatory shear of amplitude 70 and phase (j> is 
started at later time t s , so that 7(f) = @(t — t s )Re {70 exp [i(4> + wi)]}; (in) this 
is maintained up to (or beyond) a time t at which point the stress is measured. 
Using the linear constitutive equation @, we obtain 

a(t) = Re Liw /' e i{M,) G(t - t', t') d£ + l0 e i{ - Ms) G(t - t s , t s ) 

Jt a 

This unambiguously defines a time-varying viscoelastic spectrumF] as 



G*(lu, t, t s ) = iu> 



f* e'^-^Git - t', t') dt' + e-^-^Git - t s , t s ) (12) 



A similar compliance spectrum, J*(u,t,t B ) can be defined by exchanging stress 
and strain in this protocol. 

Since it depends on two time arguments as well as frequency, G*(u,t,t B ) is a 
somewhat cumbersome object. However, simplifications can be hoped for in the 
limit oj(t — t s ) S> 1. In the TTI case, this condition eliminates simple transients, 
and allows one to relate G*(u) to the Fourier transform of G(t) (see (§)). Corre- 
sponding simplifications are certainly not guaranteed in the absence of TTI. How- 
ever, the transient dependence on t s may become negligible^] when u>(t — t s ) S> 1, 
in which case we have 

G*(u,t,t B )^G*(uj,t) (13) 

giving a viscoelastic modulus that depends only on the measurement time t. If, 
in addition, the time evolution of the underlying material properties is negligible 
on the timescale of one oscillation, then G*(u>,t) may obey the relation 

poo 

G*{lj, t)= iu e- wt G(t', t) dt' (14) 







Similarly, for u>(t — t s ) ^> 1 the compliance spectrum may become t s -independent, 
J*(u,t,t s ) — > J*(u,t), and may be related to the step stress response via 



00 

iujt 



J*(u,t)=iu e~ lu)t J(t',t)dt' (15) 





Finally, G* (u, t) and J* (to, t) may obey the conventional reciprocal relation G* (u, t) 
= 1/ J*(u,t). Indeed, we shall find that all the above simplifying relationships 
are true for the SGR model studied below. As discussed at the end of Sec. ||, they 
may also hold more generally for systems with what we term there "weak long 

4 Note that in principle, to identify by experiment the real and imaginary parts of G* for a 
particular ui, t, t s one would require the experiment to be repeated for two different phases 4> (e.g. 
pure sine and cosine deformations). A more common procedure is, of course, to maintain the 
oscillatory strain over many cycles and record the "steady state" amplitude and phase response 
of the stress. For systems without TTI the latter are not uniquely defined. Only when material 
properties vary slowly enough will this give a definite result; whenever it does, it will coincide 
with The required conditions are considered, for the SGR model, below. 



'For the SGR model, an additional requirement is that u;t B 3> 1; see Sec. 3.1.2 below 
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term memory" . However, we do not have a rigorous proof for this. Pending such 
a proof, the above simplifications remain hypotheses needing explicit verification 
for any constitutive model. Experimenters should likewise beware that, for sys- 
tems without TTI, such textbook relationships between the oscillatory and step 
strain response functions cannot be assumed, but should be empirically verified, 
for each system studied. This prima facie breakdown of conventional linear vis- 
coelastic relationships in ageing systems was emphasized by [Struik (197§| ) , though 
he argued that they are recovered in sufficiently 'short-time' measurements. It 
does not (as Struik seems to suggest) extend necessarily to breakdown of linear 
superposition itself, which survives in the form of (S). 

2.8 Steady State Response: The Flow Curve 

Consider now the ultimate state of a material with TTI long after an infinitesimal 
step stress of amplitude an has been applied. The ultimate deformation may 
involve a limiting strain 7 = croJ(t — > 00), in which case the steady state (linear) 
Hookean elastic modulus is Gqq = 00/7. Alternatively, the ultimate state may 
involve a limiting strain rate, in which case the zero-shear viscosity is r/ = a/j. 
However, neither outcome need occur. If, for example, one has "power law creep", 
i.e., J(t) ~ t y with < y < 1, the material has both zero modulus (infinite 
compliance) and infinite viscosity in steady state. There is no rule against this, 
although it does require nonanalyticity of G*(uj) at small frequencies ( [Sollich 
et al., 1997|; Sollich, 199Sj ), such that T mfK is infinite. 



What if the stress amplitude is larger than infinitesimal? The ultimate steady 
state can again be that of a solid, a liquid, or something in between. In cases 
where a liquid- like response is recovered, it is conventional to measure the "flow 
curve" , which is a steady state relationship between stress and strain rate: 

0"ss = 0(7) (16) 

In many materials, the following limit, called the yield stress 

0(7 - 0) = a y (17) 

is nonzero|] Note, however, that the presence of nonzero yield stress does not 
necessarily imply a finite Hookean modulus Goo', for a < a y , the material could 
creep forever, but at an ever decreasing rate.[] Nor does the absence of a finite 
yield stress imply a finite viscosity; a counterexample is the power law fluid, 
for which a ~ j p . This has a y = but, for p < 1, infinite viscosity n = 

lim^_>oer(7)/7- 

We now turn to materials without TTI. For these, no meaningful definition of 
"steady state response" exists in general. However, in the SGR model considered 
below, TTI is restored for nonzero 7 flSollich et al., 19971 Sollich, 199SQ , and this 



6 T he experimental ex istence of a true yield stress, as denned by this limit, is debat- 
able ( Barnes et al., 1989| ); behaviour closely approximating this is, however, o ften reported . 



Note that our definition of yield stress, from the flow curve, is unrelated to that of Struik (197g| ) 



who defines a 'tensile yield stress' at constant strain rate. 

Alternatively, it could reach a steady strain 7 that is not linear in a even as a — » 0. 
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may be generic for certain types of ageing QSollich et al., 19971; Sollich, 1998| 



Bouchaud and Dean, 1995|; Kurchan, 1999|) . If so the flow curve, including the 
value of the yield stress a y (but not the behaviour for a < a y ) remains well-defined 
as a steady-state property. 

3 Ageing 

So far, we have set up a general framework for describing the rheological proper- 
ties of systems without TTI. Time translation invariance can be broken, in a triv- 
ial sense, by the transients that any system exhibits during equilibration. We now 
consider how such transients can be distinguished from ageing proper. To focus 
the discussion, we consider the linear step strain response function G(t — t w ,t w ). 
The other response functions introduced above can be treated similarly. We de- 
fine ageing (of the step strain response) as the property that a significant part 
of the stress relaxation takes place on timescales that grow with the age t w of the 
system. If ageing is present, then in order to see the full stress relaxation we need 
to allow the time t at which we observe the stress to be much larger than the 
time t w at which the step strain has been applied. Formally, we need to consider 

lim G(t-t w ,t w ) (18) 

t— *oo 

at fixed t w . On the other hand, if there is no ageing, then the full stress relaxation 
is "visible" on finite timescales. This means that as long as At = t — t w is large 
enough, we observe the full stress relaxation whatever the age t w of the system at 
the time when the strain was applied. Formally, we can take t w to infinity first, 
and then make At large, which amounts to considering 

lim lim G(At,t w ). (19) 

At— >oo t„—>oo 

In the absence of ageing, the two ways ( |l8| ) and dl9| ) of measuring the final extent 
of stress relaxation are equivalent, and we have 

lim G(t - t w , t w ) = lim lim G(At,t w ). (20) 

f^oo At— >oof w ^oo 

If the system ages, on the other hand, this equality will not hold: the right-hand 
side allows only for the decay of stress by relaxation modes whose timescale does 
not diverge with the age of the system, and thus attains a limit which includes 
elastic contributions from all modes that do have age-related timescales. It will be 
different from the left-hand side, which allows for relaxation processes occurring 
on all timescales, and thus attains a limit in which only completely non-decaying 
modes contribute. We therefore adopt the definition that a systems ages if at least 
one of its response functions violates (|20|). By contrast, we refer to deviations 
from TTI in other systems (for which all significant relaxation processes can 
essentially be observed on finite timescales) as transients. We discuss this point 



further in the context of the SGR model in Sec. 6.1.1 



Systems that violate (|20|) are referred to as having "long term memory" (puglian 



dolo and Kurchan, 1995; Bouchaud et al., 199S; Cugliandolo and Kurchan, 1993 ). 



11 



They can be further subdivided according to the strength of this memory. To 
illustrate this distinction, imagine applying a (small) step strain to a system at 
time to and switching it off again at some later time t\. The corresponding stress 
at time t > ti is proportional to G(t — to, to) — G{t — t\, t\). If this decays to zero 
at large times t, that is, if 



lim [G(t - t , t ) - G(t - ti,h)] = (21) 



[and ( |20|) is violated] then we say that the system has "weak long term mem- 
ory", otherwise it has "strong long term memory".^] Although the weakness 
condition (21) does not hold for all response functions in all ageing systems, it 
seems rather natural to expect it, in the rheological context, for most materials 
of interest. Indeed, a system with weak long term memory eventually forgets 
any perturbation that was only applied to it during a finite period. Thus, the 
treatment of a sample directly after it has been prepared (by loading it into the 
rheometer, preshearing, etc.) will not have a strong impact on the age-dependence 
of its rheological properties. This is the usual experience, and is obviously needed 
for the reproducibility of experiments results; likewise, it means that one can hope 
to make theoretical predictions which are not sensitive to minute details of the 
sample preparation. For the SGR model, any long term memory is indeed weak 



(as shown in Sec. 6.1.1 below); we consider this an attractive feature. Note in any 
case that a rheological theory for systems with strong long term memory might 
look very different from the SGR model. 

We have defined ageing as the property that a significant part of the stress 
relaxation G(t — t w ,t w ) takes place on timescales that grow with the age t w 
of the system. In the simplest case, there is only one such growing timescale, 
proportional to the age of the system itself. The (ageing part of the) stress 
relaxation then becomes a function of the scaled time difference (t — t w )/t w . 
We will encounter such simple ageing behaviour in the glass phase of the SGR 
model, which is discussed below. More complicated ageing scenarios are possible, 
however: There may be several timescales that grow differently with the age of 
the system. This can be represented as 

G(t- ) = £&[fti(*)/M*w)] (22) 

i 

where the functions hi{t) define the different diverging timescales. If there is 
only a single term in the sum, with h{t) = t, then the simplest ageing scenario 
(shown by the SGR model) is recovered. On the other hand, for h(t) = exp(t/ro) 
(where to is a microscopic time) one has TTI. The more general form h(t) = 
exp[(t/ro) 1_/1 ] interpolates between these two limiting cases (and, for i — i w <§C t w , 
gives Struik's general 'time waiting-time superposition principle' ( jBtruik, 1978| )). 

8 There is a slight subtlety with the definition of long term memory for the linear step stress 
response. Eq. (blj) , applied literally to J(t — t w ,£w), suggests that even a Newtonian fluid with 
J(t — t w ,iw) ~ t — t w has strong long term memory, because its strain "remembers" stress 
applications in the arbitrarily distant past. This is clearly undesirable as a definition. The 
problem can be cured by "regularizing" the step stress response: one simply considers the 
material in question "in parallel" with an elastic spring with infinitesimal modulus. 
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Otherwise, Cugliandolo and Kurchan have shown under fairly mild assumptions 
that (22) is the most general representation of the asymptotic behaviour of 
step response and correlation functions in systems with weak long term mem- 
ory flCugliandolo and Kurchan, 1994 ). 

Let us return now to the status of Eqs. (13,14,1^). (These concern the lack 
of ^-dependence in G*(u,t,t s ), the Fourier relationship between frequency and 
real-time spectra, and the reciprocity between G* and J*.) As stated in Sec. 2.7 
these equations have no general validity for systems without TTI. Indeed, one can 
easily construct theoretical model systems with strong long term memory which 
violate them. On the other hand, we speculate that systems with weak long term 
memory will generically have the properties (|lj| 14,15|). Plausibility arguments 
can be given to support this hypothesis ( [Fielding, 200C| ), but these do not yet 
amount to a proof. The cautionary remarks at the end of Sec. 2/7 therefore still 
apply. 



4 The SGR model 

The SGR model is a phenomenological model which captures many of the ob- 
served rheological properties of soft metastable materials, such as foams, emul- 



sions, slurries and past 


es flMackley et al., 1994; Ketz et al., 1988; Khan et al., 


1988; Mason et al., 19951 Panizza et al., 1996|; Hoffmann and Rauscher, 1993; 


Mason and Weitz, 1995 


). It is based upon Bouchaud's trap model of glassy dy- 



namics, with the addition of strain degrees of freedom, and the replacement of the 
thermodynamic temperature by an effective (noise) temperature. It incorporates 
only those characteristics deemed common to all soft glassy materials (SGM's), 
namely structural disorder and metastability. We now review its essential fea- 
tures. 

We conceptually divide a macroscopic sample of SGM into many mesoscopic 
elements. By mesoscopic we mean large enough such that the continuum variables 
of strain and stress still apply for deformations on the elemental scale, and small 
enough such that any macroscopic sample contains enough elements to allow the 
computation of meaningful "averages over elements". We then assign to each 
element a local strain I, and corresponding stress kl, which describe deformation 
away from some local position of unstressed equilibrium relative to neighbouring 
elements. The macroscopic stress of the sample as a whole is defined to be {kl), 
where () denotes averaging over elements. Note that, for simplicity, (shear-) 
stress and strain are treated as scalar properties. The model therefore does not 
predict, or allow for, the various normal stresses which can arise in real materials 



undergoing nonlinear shear ( Doi and Edwards, 1986 ). 



For a newly prepared, undeformed sample, we make the simplest assumption 
that I = for each element. (Physically, of course, (I) = would be sufficient and 
is indeed more plausible.) The subsequent application of a macroscopic strain at 
rate j causes each element to strain relative to its local equilibrium state and ac- 
quire a non-zero I. For a given element, this continues up to some maximal strain 
l y , at which point that element yields, and rearranges into a new configuration of 
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local equilibrium with local strain I = 0.^] Under continued macroscopic strain- 
ing, the yielded element now strains relative to its new equilibrium, until it yields 
again; its local strain (and stress) therefore exhibits a saw-tooth dependence upon 
time. 

The simplest assumption to make for the behaviour between yields is that 
7 = I: the material deformation is locally affine (Doi and Edwards, 1986). Yield 
events apart, therefore, the SGR model behaves as an elastic solid of spring 
constant k. Yields confer a degree of liquidity by providing a mechanism of stress 
relaxation. 

Although above we introduced yielding as a purely strain-induced phenomenon, 
we in fact model it as an "activated" process ( Sollich et al., 1997; Sollich, 1998 ). 
We assume that an element of yield energy E = ^kl 2 , strained by an amount 
I, yields with a certain rate; this defines the probability for yielding in a unit 
time interval. We write this rate as r _1 , where the characteristic yield time 
r = roexp (E — ^kl 2 )/x is taken to be the product of an attempt time and 
an activation factor which is thermal in form. This captures the strain-induced 
processes described above since any element strained beyond its yield point will 
yield exponentially quickly; but it also allows even totally unstrained elements 
to yield by a process of activation over the energy barrier E. These activation 
events mimic, within our simplified model, nonlinear couplings to other elements 
(the barrier heights depend on the surroundings, which are altered by yield events 
elsewhere). A more complete model would treat these couplings explicitly. How- 
ever, in the SGR model, which does not, x is regarded as an effective "noise" 
temperature to model the process. Because the energy barriers are (for typical 
foams, emulsions, etc.) large compared to the thermal energy ksT, so are the 
energy changes caused by these nonlinear couplings, and so to mimic these, one 
expects to need x of order the mean barrier height Note that the SGR 

model treats "noise-induced" yield events (where the strain is much below the 
yield strain l y , i.e., where \kl 2 <C E) and "strain-induced" yield events (where 
\kl 2 ~ E) in a unified fashion. We will nevertheless find it useful below to 
distinguish between these two classes occasionally. 

The disorder inherent to SGM's is captured by assuming that each element 
of a macroscopic sample has a different yield energy: a freshly yielded element is 
assigned a new yield energy selected at random from a "prior" distribution p(E). 
This suggests the following alternative view of the dynamics of the SGR model, 
which is represented graphically in Fig. [2|. Each material element of a SGM can 
be likened to a particle moving in a landscape of quadratic potential wells or 
"traps" of depth E. The depths of different traps are uncorrelated with each 



9 This ignores possible "frustration" effects: an element may not be able to relax to a fully 
unstrained equilibrium position due to interactions with neighbouring elements. Such effects 
can b e incorporate d into the model, but are not expected to affect the results in a qualitative 
way flSollich, 1998| ). 

10 Whether i t is fully consistent to have a nois e temperature x 3> fesT is a debatable feature 
of the model (Sollich et al., f997; Sollich, f998); however, we think the results are sufficiently 
interesting to justify careful study of the model despite any uncertainty over its interpretation. 
It is also intriguing to note that similar "macroscopic" effective temperatures (which remain 
nonzero even for ksT — » 0) have recently been found in other theo ries of out-of-equilibrium 
systems with slow dynamics (Kurchan, 1999; Cugliandolo et al., 1997). 
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other and distributed according to The bottom of each trap corresponds 

to the unstrained state I = 0; in straining an element by an amount I, we then 
effectively drag its representative particle a distance \kl 2 up the sides of the 
trap, and reduce the effective yield barrier height (E — > E — ^kl 2 ). Once the 
particle has got sufficiently close to the top of its trap (E — \kl 2 w x), it can hop 
by activated dynamics to the bottom of another one. This process corresponds 
to the yielding of the associated material element. In the following, we shall 
use the terminology of both the "element picture" and the "particle picture" as 
appropriate. Thus, we will refer to r = rrjexp (E — \kl 2 )jx as either the yield 

or relaxation time of an element, or as the lifetime of a particle in a trapQ The 
inverse of r is the rate at which an element yields/relaxes or a particle hops. 
However, we normally reserve the term yield rate or hopping rate for the average 
of these rates over the whole system, i.e., over all elements or particles. This 
quantity is denoted Y and will occur frequently below. 

A specific choice of p(E) is now made: p(E) = (l/x g )exp(—E/x g ), where 
x g = (E) is the mean height of a barrier chosen from the prior distribution p(E). 
As shown by Bouchaud (1992| ), the exponential distribution, combined with the 



assumed thermal form for the activated hopping, is sufficient)^] to give a glass 
transition in the model. The transition is at x = x g and divides the glass phase 
[x < x g ), in which weak ergodicity breaking occurs, from a more normal phase 
(x > x g ). In the glass phase, the Boltzmann distribution (which is the only 
possible steady state for activated hopping dynamics, in the absence of strain), 

PeqOE) oc p(E) exp(E/x) (23) 

is not normalizable: thus there is no steady state, and the system must age with 
time. (The converse applies for x > x g : there is then a unique equilibrium state, 
which is approached at long times. Hence ageing does not occur, though there 
may be transients in the approach to equilibrium.) Apart from our use of an 
effective temperature x, the only modification to Bouchaud's original model of 
glasses lies in our introduction of dynamics within traps coupled to strain. 

It may appear suspicious that, to obtain a glass transition at all, an expo- 



nential form of p(E) is required ( Bouchaud, 1992 ). In reality, however, the glass 



transition is certainly a collective phenomenon: the remarkable achievement of 
Bouchaud's model is to represent this transition within what is, essentially, a 
single-particle description. Thus the chosen "activated" form for the particle 
hopping rates, and the exponential form of the trap depth distribution, should 
not be seen as two independent (and doubtful) physical assumptions, but viewed 
jointly as a tactic that allows glassy dynamics to be modelled in the simplest 
possible way flSollich et al., 1997|; Sollich, 1998|) . 



From now on, without loss of generality, we choose units so that both x g and 
k are unity. This means that the strain variable I is defined in such a way that 



11 Because of this lack of correlation, it does not make sense to think of a particular spatial 
arrangements of the traps. 

12 Sometimes this will be further abbreviated to "lifetime of a trap" or "lifetime of an element" . 

13 This is sufficient, but it is necessary only that the given exponential form be approached at 
large E. 
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an element, drawn at random from the prior distribution, will yield at strains of 
order one. Since the actual value of the strain variable can be rescaled within 
the model (the difference being absorbed in a shift of k), this is purely a matter 
of convention. But our choice should be borne in mind when interpreting our 
results for nonlinear strains, given below: where strains "of order unity" arise, 
these are in fact of order some yield strain L, which the model does not specify, 
but which may in reality be a few percent or less. In addition we choose by 
convention To = 1; the timescale in the SGR model is scaled by the mesoscopic 
"attempt time" for the activated dynamics. The low frequency limit, which is 
the main regime of interest, is then defined by ujtq = wCl, Note that, without 
our choice of units, (E) = 1 so that we expect the interesting physics to involve 
x ~ 1. 



E-l/2kV 




Figure 2: Dynamics of the SGR model. A representative particle (1) may hop 
out of its trap by activated hopping (1 — ► 2). It enters the new trap in a state 
of zero local strain (I = 0); application of strain A7 raises its energy (2 — ► 3) 
making a subsequent hop (3 — > 4) more likely. Note that the relative horizontal 
displacement of the quadratic potential wells (traps) is arbitrary; each has its 
own independent zero for the scale of the local strain I. 



4.1 Constitutive Equation 

The SGR model is exactly solved by two coupled constitutive equations flSolfich, 
1998| ), the first of which expresses strain as an integral over stress history, while 
the second embodies the conservation of probability. We assume that the sample 
is prepared (in a known initial state of zero stress and strain) at time zero and 
that a time dependent macroscopic strain 7(f) is applied thereafter, so 7(2) = 
for t < 0. The constitutive equations are then 

a(t) = 7 (t)G (Z(t, 0)) + f [ 7 (t) - 7 (0] Y(t')G p (Z(t, t'))dt> (24) 

J 
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1 = G (Z(t, 0)) + f Y(t')G p (Z(t, t'))dt' (25) 
J o 

In these equations 

Z(t,t') = £ exp ( [ 7 (f") - 7 (0] 2 /2s) di" (26) 
and G P (Z) and Gq{Z) obey 

G P (Z) = p(£) exp (-Ze~ E/x ) dE (27) 

G (Z)=j P {E)exp(-Ze- E/x )dE (28) 

where Po(E) is the probability distribution for the yield energies (or trap depths) 
in the initial state of preparation of the sample at time t = 0. We return below 
(Sec. 5.1) to the issue of how to choose this initial state. 

These equations can be understood by viewing yielding as a "birth and death" 
process: each time an element yields it dies and is reborn with zero stress, and 
with a yield energy selected randomly from the prior distribution p(E). The (av- 
erage) yield rate rate at time t' is Y(t'); the birth rate at time t' of elements of 
yield energy E is therefore Y(t')p(E). The proportion of these which survive with- 
out yielding until time t is found as exp [—Z{t,t') /t(E)\ where t(E) = exp(E/x) 
is the (mean) lifetime that an unstrained element with yield energy E would 
have. The expression ( [26] ) for Z(t, t') reflects the fact that an element that 
has last yielded at time t' and has a yield energy E will have a yield rate of 
t(E)^ 1 exp f [ 7 (t") — 7 (i')] 2 /2z) a t time t". Here the exponential factor accounts 
for the lowering of the yield barrier by strain applied since the element last yielded 
(see Fig. |2|) . Note that this factor is unity under conditions where the local strain 
is everywhere negligible, in which case Z(t, t') — ► t — t' (we return to this point 
below). More generally, Z(t,t') can be thought of as an effective time interval 
measured on an "internal clock" within an element, which allows for the effect of 
local strain on its yield rate, by speeding up the clock. This speeding up effect, 
which describes strain-induced yielding, is the only source of nonlinearity within 
the SGR model. 

According to the above arguments, the number of elements of yield energy E, 
present at time t, which were last reborn at time t' is 

P(E, t, = Y{t')p{E) exp [-Z(t, t')/r(E)} (29) 

Such elements each carry a local strain 7 (i) — j(t') and so the net contribution 
they make to the stress at time t is 

s(E, t, = [ 7 (t) - 7 (0] Y(t')p(E) exp [-Z(t, t')/r{E)] (30) 

Integrating these expressions over t' from to t and adding terms representing 
the contribution from elements which have survived from t = without yielding 
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at all, we get respectively the number P(E, t)dE of elements at time t with yield 
energies between E and E + dE: 



P{E,t) = P (E)exp \-Z(t,0)e- E/x 



t 

/ P(E,t,t')dt' 
'o 



(31) 



and the corresponding stress contribution s(E, t)dE at time t from such elements: 



s(E,t) =~f{t)P {E)exp \-Z(t,0)e- E/x ] + / s(E,t,t')dt' 



(32) 



Integrating (|3l| ) and (|32| ) over all yield energies E, we finally recover our consti- 
tutive equations (p4|) and (p5|) respectively. Below we will return to these two 
quantities, which will sometimes be expressed instead as a function of the life- 
time t(E) = exp(E/x) of an unstrained element with yield energy E, so that 
P(r,t)dT = P(E,t)dE and likewise for s. Note that, because E > 0, these distri- 
butions are nonzero only for r > 1. We will not write this restriction explicitly 
below. 

Finally, the following alternative form of the first constitutive equation ( p4j ) 
is sometimes useful: 



a(t)= 7 (t)- / 7(W)G^(M'))^ 
Jo 



This is obtained by substituting 
Z(t,t') is again replaced by i — i'. 



into <\2 



(33) 

In the limit of small strains, 



4.2 Rheological Properties of the SGR Model 



Solution of the constitutive equations (24, 25) is relatively straightforward un- 
der conditions where TTI applies. Here we recall the main results thereby ob- 
tained ( ftollich et al., 1997]; Sollich, 1998| ). 



4.2.1 Linear Spectra 



A regime of linear rheological response arises whenever the effects of strain on 
the effective time interval Z(t, if) is small. This requires that the local strains 
in each element remain small; in oscillatory shear, where j(t) = jqc 1 ^, this is 
satisfied at low enough strain amplitudes 70 for any finite frequency u>. (The 
same is not true in steady shear flow; we return to this in Sec. 4.2.2 below.) In 
the linear regime, the model's internal dynamics are independent of the imposed 
deformation: the elements' lifetimes are, to order 70, strain-independent. In the 
constitutive equations, Z(t, t') can then be replaced by the time interval t — t' 
(there is no strain-induced yielding). 

As described in Sec. p. 61 above, the conventional definition of the linear vis- 



coelastic spectra G' (to) , G" (lo) (Eqs. 10,11), requires not only linearity but also 
TTI. Thus they are well-defined only for an equilibrium state; in the SGR model, 
the latter exists only for x > 1. But even at x > 1 these spectra show interesting 
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power law dependencies at low frequency ^; these are summarized as follows (the 
prefactors are omitted, but discussed by Sollich et al. (1997); Sollich (1998|) ): 



G" 

a 



oc 
oc 



uj for 2 < x, 
uj 2 for 3 < x, 



oc uj x 1 for 1 < x < 2 
oc uj x ~ 1 for 1 < x < 3 



(34) 



Throughout its glass phase (x < 1) where the SGR model violates TTI, we 
must study instead the time dependent spectra G*(oj,t,t s ) as defined in Sec. 2.7 



above; this is done in Sec. 6.1 below. An alternative, explored by Sollich et al 



(1997|); Sollich (1998|); Evans et al. (1999j ) is to observe that TTI can be restored 
even for x < 1 by introducing a cutoff E max in the trap depth distribution p(E). 
This gives interesting predictions for x < 1: for example, one finds G'(oj) ~ lu 1 ~ x , 
for T- l (E max ) <C uj < 1 QSollich et al., 1997|; Sollich, 1998Q . However, the role 
of this cutoff is to bring all ageing processes to a halt after a large finite time 
of order r(£' max ); formally there is no long term memory. Since in the present 
work we want to study the ageing regime itself, we assume instead that E mavi is 
infinitely large, so that for x < 1, ageing continues indefinitely. 



4.2.2 Flow Curve 



The flow curve was defined in Sec. 2^ as the nonlinear stress response cr(j) to a 
steady strain rate 7. For the SGR model, it shows the following scalings: 



(7 
O 
(7 



OC 
OC 
OC 



(35) 



7 for x > 2 
7 1 " 1 for 1 < x < 2 
<7 y oc 7 1 ~ x for x < 1 

Here 7 <C 1 is assumed; prefactors are discussed by Sollich (1998j ). The flow 
curve exhibits two interesting features which are explored more fully in Sees. |6.2.2 
Firstly, for x 



an d f7~2~l 
(gdlicl 



< 1 there is a yield stress <J y (x) (whose value is plotted 
in ( jSollich, 1998| ) ) . A linear response regime exists at <r C d y ; ageing can occur 
for all a < a y . For a > a y the system achieves a steady state, and ageing no 
longer occurs. This is because any finite flow rate, however small, causes strain- 
induced yielding of elements even in t he deepest traps -P| Thus the agein g process 



is curtailed or "interrupted" by flow ( Sollich et al., 1997; Sollich, 1998|) ; the flow 



curve is well-defined (and independent of the choice of Po m the initial state) 
even in the glass phase. The second interesting feature is that, for 1 < x < 2 
(where ageing is absent) there is no linear response regime at all in steady shear: 
however small the applied stress, the behaviour is dominated by strain-induced 
yielding. There is an anomalous (power law) relation between stress and strain 



rate, and an infinite zero-shear viscosity (cf. Sec. 2.8 above). This also shows up 
in (34), where rj = lim^^o G"(u)/u is likewise infinite. 

14 Here and throughout this paper, "low frequency" in the SGR model means, uj 1, that is, 
frequencies small compared to the mesoscopic attempt rate for activated hopping Tq = 1 (in 
our chosen units). 

15 The time required to yield, with a steady flow present, is only power law, rather than 
exponential in E. 
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5 Ageing in the SGR model 



In this section we discuss some general features of ageing in the SGR model; in 
subsequent ones, we explore the rheological consequences of these phenomena. 

5.1 Initial Preparation of Sample 

As noted above, to solve the constitutive equations (24,25[) the initial distribution 
Pq(E) of yield energies or trap depths at time zero must be specified. Since we 
are largely interested in the rheological properties of the glass phase (x < 1), 
for which no steady-state distribution of yield energies exists in the absence of 
flow, we cannot appeal to equilibrium to fix Pq(E). Instead, this should depend 
explicitly on the way the sample was prepared. For simplicity, we choose the 
case where Po(E) = p(E); this is equivalent to suddenly "quenching" the noise 
temperature x, at time zero, from a very large value {x > 1) to a value within 
the range of interest. We refer to it as a "deep quench". 

The question of whether or not a deep quench is a good model for the sample 



preparation of a SGM remains open ( Somen et al., 1997; Sollich, 1998| ); since x 



is not truly a temperature, it is not clear exactly how one would realize such a 
quench experimentally.p| However, we expect that most interesting aspects of 
ageing behaviour are not too sensitive to the initial quench conditions Po(E), so 
that a deep quench is indeed an adequate model. A study of the effect of quench 



depth on the results for the SGR model is summarized in App. A. 4; we find 
independence of quench depth so long as the final noise parameter x is not too 
small.|^| More generally, a degree of insensitivity to the initial quench conditions 
is consistent with the weak long term memory scenario; a system whose response 
decays with a relaxation time of order its age will typically lose its memory of 
the initial state by a power law decay in time. This can then easily be swamped 
by larger, Pq independent contributions, as indeed occurs in most regimes of the 



SGR model (App.|AJ). 

Following the initial preparation step, subsequent time evolution of the rhe- 
ological response is, within the glass phase, characterized by an ageing process. 
To allow simpler comparisons with the non-ageing (but still slow) dynamics for 
1 < x < 2, below we shall also consider a similar quench from large x to values 
lying in this range. 

5.2 Ageing of the Lifetime Distribution 



We now (following [Bouchaud (1992 ) and Monthus and Bouchaud (1996| )) dis- 



cuss in detail the way ageing affects the lifetime distribution (or equivalently the 
distribution of particle hopping rates) within the SGR model. 

16 One argument in its favour is that it this choice minimizes the information content (maxi- 
mizes the entropy) of the initial distribution Po; it is therefore a legitimate default choice when 
no specific information about the preparation condition is available. 

17 More precisely, if the "deep quench" specification is altered to one in which, at time zero, 
the system is quenched from equilibrium at xq > 1 to its final noise temperature x, the leading 
results are independent of xo so long as the final x value obeys x > 1/(2 — 1/xq). Note that this 
condition is never satisfied for x < 1/2. 
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We ignore the presence of a strain; the following results apply when there 
is no flow, and in the linear response regime, where strain-induced hops can be 
ignored. Under such conditions, the hopping rate Y(t) is a strain-independent 
function of time, and is readily found from ( p5|) by Laplace transform. This is 
done in App. |A.2| . For the case of a deep quench (as defined above), the exact 
asymptotic forms of Y are as follows: 

Y(t) = ^— ^ for x > 1 

x 

Y{t) = W) fOT X = 1 (36) 

Y(t) = ** for x<\ 

XL{X)T{1 — X) 

where T(x) is the usual Gamma function. These results assume t 3> 1, which we 
will usually take to be the case from now on (since timescales of experimental 
interest are expected to be much longer than the mesoscopic attempt time tq = 1). 
Note that the late-time asymptotes given here are subject to various subdominant 
corrections (see App. |A.2| ), some of which are sensitive to the initial state of 
sample preparation.^ 

A closely related quantity to the hopping rate Y is the distribution of yield 



energies P(E,t) - which obeys (31) - or equivalently the lifetime distribution 
P(t, t). As previously pointed out, in the absence of strain, the only candidate for 
a steady state distribution of yield energies P eq (E) is the Boltzmann distribution: 
P eq (E) oc p{E) exp(£'/x), which translates to P C q(T) = P cq (E)dE/dr cx t~ x ; in 
either language, the distribution is not normalizable for x < 1, leading to broken 



TTI in the model ( [Bouchaud, 1992|) . 



Let us therefore consider a deep quench at time t = 0, and define the proba- 
bility distribution for trap lifetimes P(r, i w ) as a function of the time t w elapsed 
since sample preparation. (In Sec. [], we will identify t w with the onset of 
a step strain.) The initial lifetime distribution, P(r, 0), describes a state in 
which the trap depths are chosen from the prior distribution P(E,0) = p(E); 
just after a quench to temperature x the distribution of lifetimes is therefore 
P(t, 0) x p(E)dr I 'dE x t~( 1+x \ Thereafter, by changing variable from E to r 
in (|3l|) , we find the following approximate expressions for P(t, i w ) 

P(r,i w ) ~ xY(t w )rp{T) for r < t w and t w » 1 , . 

P(t, t w ) ~ xY(t w )t w p{r) for r»t w and t w »l [ ) 

For a quench temperature above the glass point (x > 1), P(r,t w ) exhibits a 
transient decay; as t w — > oo, we find (using the results in fl36"D) that P(r,t) — > 
-f'eq(' r ) = (1 — ^)t _3: , as expected. The nature of the approach to the long time 
limit is illustrated schematically in Fig. |3|(a); the final distribution has most of its 

18 For a quench from initial noise temperature xo, the relative order of the affected subdominant 
terms becomes t~ x ^ 1 ~ 1 ^ x °\ Thus, unless one quenches from a point that is itself only just above 
the glass transition, or to a point that has x only just above zero, the exact specification of the 
initial state is unimportant at late times. 
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weight at r = 0(1), consistent with the fact that the hopping rate ( |36[ ) is itself 
O(l) in this phase of the model. 
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Figure 3: Schematic evolution of the relaxation time distribution (a) above the 
glass transition; (b) below it. The first shows a transient decay onto a steady 
state, the second shows ageing behaviour. The curves lie in order of increasing 
t w at the bottom of each figure. 

For x < 1, in contrast, P(r,t w ) evolves as in Fig. ||(b); the limit of P(r, i w ) 
is zero for any finite r as t w — » oo. Hence, the proportion of elements having 
yield time of order unity tends to zero as i w — > oo; the bulk of the distribution's 
weight is at r ~ t w . 19 This is consistent with the decay of the hopping rate as 
a power law of t w , and with the idea that, in a system undergoing ageing, the 
characteristic relaxation time is typically of the order the age of the system itself. 



5.3 Higher Moments 

The above analysis focuses on the time-evolution of the distribution of elements' 
lifetimes, which is the usual quantity of interest in the formal analysis of ageing 
effects ( Bouchaud and Dean, 1995; Bouchaud et al., 1998 ). Indeed, the latter 
are usually attributed to the divergence of the normalization integral, or zeroth 
moment, of feq( T ) (undefined, within the Boltzmann distribution, when x < 1). 
Formally, however, one can consider a series of critical x values, x n = n + 1, below 
each of which the nth moment of P eq becomes undefined ( [Evans et al., 199 9 



19 | More formally, for x < 1, we have lim^^oo P(t, t w )d,T = for any b > 1, while for any 
a < 1 < b we have instead lim^^oo J^" P(r, t w )dr = O(l). 
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Odagaki, 1995). For n > this does not lead to ageing, in the sense defined in 
Sec. H| above, but can lead to anomalous, slow time evolution in any experimental 
measurement that probes the nth moment. For example, in Sec. |6. 1 .3 below, we 



discuss the time-evolution of the distribution of stresses borne by elements in a 
steady-shear startup experiment. In steady state, the stress carried by an element 
whose lifetime is r is of order jr. If P(t) = P e q(T) and is unperturbed by flow 
(as a linear response analysis would assume), then the zero-shear viscosity is of 



order J tP ccl {t)(1t, which diverges for x < 2 (see Sec. ^8 above 



6 Rheological Ageing: Imposed Strain 

In this and the next sections, we describe our new rheological results for the SGR 
model. We focus particularly on rheological ageing, which occurs in the glass 
phase (x < 1); however, several new results for 1 < x < 2, including anomalous 
transient behaviour, are also presented. The case x = 1, which divides these 
regimes, shows its own especially weak (logarithmic) form of ageing and is, where 
necessary, treated separately below. 

For simplicity, we consider (for all x values) only the idealized route to sample 
preparation described in Sec. above: the system is prepared at time t = 
by means of a deep quench, so that Go(Z(t, 0)) = G p (Z(t, 0)) in the constitutive 



equations (24, 25). Note that these constitutive equations for the SGR model are 
more readily solved to find the stress response to an imposed strain, rather than 
vice-versa. Accordingly, we focus first on strain-controlled experiments and defer 
to Sec. our analysis of the stress-controlled case. 

6.1 Linear Response 



As described in Sec. 4.2.1 above, when local strains are negligible, the SGR 



model displays a linear response regime. The effective time interval Z{t, t') in 



Eqs. (24,25) becomes the actual time interval t — t', and the hopping rate Y(t') a 
strain-independent function of time. For the deep quench considered here, Y(t') 
assumes the asymptotic forms summarized in (p6|). The stress response to any 



strain history then follows simply from (p4[), by integration. 
6.1.1 Step Strain 

For a step strain, the amplitude 70 gives the maximum local strain experienced by 
any element. The condition for linearity in this case is therefore simply 70 <C 1. 
The linearized step strain response was defined in @. It is found for the SGR 
modeQ using (|33|): 

G(t - t w , t w ) = l- C Y(t')G p (t - t')dt' (38) 



Note that by construction of the SGR model, the linear step strain res ponse is actually 



identical to the correlation fun ction defined by Bouchaud for his trap model (Bouchaud, 1992 



Monthus and Bouchaud, 199E) 
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As outlined in App. A. 3, analytic limiting forms for G(t — t w ,t w ) can be found 
when experimental timescales are large on the scale of the mesoscopic attempt 
To = 1, so that t — t w 3> 1 and t w 3> 1. In this limit we identify two distinct 
regimes: a short time interval regime t — t w <C t w an d long time interval regime 
t — t w S> t w (where the measure of "short" and "long" is not now To but t w 
itself). The limiting forms in each case depend on the value of x; our results are 
summarized in table |j. 





G(t 




) for t — t w <C t w 


G(t — t w , t w ) for t — t w » t w 


X > 1 




r(x) 


(t - t w f- x 


{X mX \t-t w f 


X = 1 




i - 


In (t - t w ) 


tw 1 




In (t w ) 


t - tw In (t w ) 


X < 1 


1 - 


i 


(t-t w y-* 


1 / tw V 


T(2-x 


)r(x) \ t w J 


r(i + x)r(i-x) \t-t w J 



Table 1: Stress response to step strain at short and long times (t-t w > 1, t > 1 
assumed). T(x) denotes the usual Gamma function. 



The asymptotic scalings apparent in the various entries of table |l| can be 
physically motivated by the following simple arguments. Upon the application of 
the strain at time t w the local strain of each element exactly follows the macro- 
scopic one, and the instantaneous response is elastic^]]: G(0, t w ) = 1. In the time 
following t w , elements progressively yield and reset their local stresses I back to 
zero. The stress remaining at t will be that fraction of elements which has sur- 
vived from t w without yielding, and hence roughly that fraction f™ t P(r, t w ) dr 
which, at time t w , had time constants greater than t — t w . Hence in measuring 
the linear response to a step strain we are probing the properties of the system 
as they were at the time of strain application. 

Using the approximate expressions given in (37) above, we have P(r, t w ) oc 
t~ x for t <§[ t w and P(r, t w ) esc t w r~( 1+:E ) for r 3> t w - This gives, for short time 
intervals (t — t w <^t w ) 

^ »n w.^ (i-iw) 1 -*-! 



/t — tw 
P(t, t w ) rfr ~ 1 - x 



tw 3? 



21 This is a general characteristic of the SGR model: whenever the macroscopic strain changes 
discontinuously by an amount Ay, the stress a also increases by A7. 
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Figure 4: Approximate curves for G(t — t w , t w ) generated using the interpolating 
formulae ( j39|]4"o| ) (dashed lines), compared to numerical data for this quantity 
(solid lines). Solid curves downwards show numerical data for G vs t — t w at 
noise temperatures x = 0.7, x = 1.0 and x = 1.3 respectively. The waiting time 
t w is 10 5 . 



and, for long time intervals (t — t w 2> t v 



P{r,t w )dT ~ 



(i - x)t w (t - t w y 



In fact these estimates approximate the numerical data already quite well. Even 
better agreement is obtained by adjusting the prefactors to fit the asymptotic 
results in table |]: 

iXxXi-tw) 1 -* 



G(t — t w ,t w ) 
G(t — t w , t w ) 



1 



1 



l-x 



T 2 {x)T(2-x)t 
{x-l)T(x)t w (t-t w y 



1 - T{x)T(x + l)r(2 - x)ti' 



l-x 



for t - t w < t w (39) 
for t - t w > t w (40) 



In the relevant time regimes, these formulae agree well with our numerical results 
(see Fig. |4|), at least over the noise temperature range to 2; they could therefore 
be used in a standard curve fitter for comparison with experimental data. In 
the limit t — t w — > oo and t w — > oo, they reproduce (by construction) the results 
shown in table [l] for x > 1 and x < 1. The logarithmic terms at the glass point 
(x = 1) can also be recovered by taking the limit x — > 1 first. 

Using the forms for G(t — t w , t w ) as summarized in table [l], and substituting 



these in Eqs. (2C,21), we see that the SGR model has short term memory for 
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(HX, (t-tj/t w 

Figure 5: Left column: age-dependent stress relaxation modulus G{t — t w ,t w ) 
against scaled time interval (t — t w )/t w (for x < 1) and time interval t — t w (for 
x > 1). Right column: 1 — G(t — t w ,t w ), plotted similarly. Shown are data for 
waiting times t w = 10 2 , 10 3 . . . 10 6 (left to right for top four graphs, right to left 
for bottom four graphs). Transients are visible in the top left figure, as follows: 
The curves coincide at short time intervals t — t w <C t w . At large t w , this regime 
accounts for more and more of the decay of G; the remaining i w -dependence is 
only through an unimportant tail. For t w — > oo, the "short time" regime extends 
to all finite values of t — t w ; one recovers the equilibrium response (shown as the 
dotted line) which decays to zero on a t w -independent timescale. Ageing is visible 
at bottom left, where the major part of the decay of G occurs on a timescale of 
order t w itself, with unimportant corrections to this scaling at early times. 
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x > 1 and weak long term memory for x < 1. Thus we expect transients for 
x > 1 and ageing for x < 1. As elaborated in Fig. || this is indeed what we 
find. More generally, these step strain results for the SGR model show some 
interesting features of rheological ageing. Consider first the behaviour above the 
glass transition (x > 1). Here the stress decay at short time intervals (t—t w <C t w ) 
depends only upon the time interval between stress imposition and measurement 
itself (t — t w ), and not on the sample age t w . This is because the traps which 
contribute to stress decay during this interval are mainly those with lifetimes 
r < t — t w ; and the population of these traps has already reached Boltzmann 
equilibrium before the strain is switched on (see Fig. ||(a)). Taking the limit 
i w — > oo at constant t — i w , (i.e., letting the system fully equilibrate before we 
apply the strain) , we recover a TTI stress relaxation function which decays to zero 
on timescales of order one (the mesoscopic attempt time). On the other hand, 
for any finite waiting time t w , the stress decay at long enough times (t — t w » t w ) 
violates TTI, since it is controlled by decay out of deep traps (r 3> t w ) which 
had not already equilibrated before t w . Note that even though this feature of 
the stress relaxation depends explicitly on t w , it is not an ageing effect according 
to our definition in Sec. I This is because the deviations from TTI and the 
dependence on t w manifest themselves at ever smaller values of G as t w becomes 
large. Equivalently, if wg ctssuniG that G(t — t w , t w ) can be measured reliably 
only as long as it remains greater than some specified value (a small fraction 
e of its initial value G(0,t w ) = 1, for example), then the results will become 
i w -independent for sufficiently large t w . 

Below the glass point (x < 1) we see true ageing, rather than anomalous 
transient effects: A significant part of the stress relaxation G{t — t w ,t w ) now 
takes place on timescales that increase with the sample age t w itself. In fact, 
in the case of the SGR model, this applies to the complete stress relaxation, 
and t w itself sets the relevant timescale: for x < 1, G depends on time only 
through the ratio (t — i w )Aw-H R is still true that stress decay during the 
interval t — t w is dominated by traps for which r < t — t w , but no longer true 
that these traps have reached Boltzmann equilibrium by time t w : in an ageing 
system such equilibrium is never attained, even for a subset of shallow traps (see 
Fig. H|(b)). Instead, the population of such traps will gradually deplete with age, 
as the system explores ever-deeper features in the energy landscape. Decay from 
these deep traps becomes ever slower; the limit t w — > oo (for any finite t — t w ) 
gives completely arrested behaviour in which all dynamics has ceased, and the 
system approaches a state of perfect elasticity (G = 1). Even in an experiment 
that can only resolve values of G above a threshold e (see above), we would detect 
that the stress relaxation becomes slower and slower as t w increases. 

The fact that G depends on time only through the ratio (t — t w )/t w is a simple 
example of Struik's 'time ageing-time superposition' principle ( [Struik, 1978 ): the 

22 This is typical, but not automatic for ageing systems; the case x = 1, for example, does not 
have it. In general, the timescale for ageing can be any monotonically increasing and unbounded 
function of £ w - There can also be parts of the stress relaxation which still obey TTI. An example 
is G(t — t w , £ w ) = gi (t — t w ) + g2 ((t — t w ) /tw) , which exhibits ageing when <?2 is nonzero , but also 
has a TTI short time part described by g\. Superpositions of relaxations with different ageing 
timescales are also possible; compare Eq. fel). 
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relaxation curves for different t w can be superposed by a rescaling of the time 
interval t — t w by the sample age. However, as mentioned previously, Struik's 
discussion allows a more general form in which the scale factor varies as t!^, with 
fi < 1. The case fi = l, exemplified by the SGR model, is the only one in which, 
even at very long times, the relaxation time does not become short compared to 
system age. 



6.1.2 Oscillatory Strain 

In an oscillatory strain, the maximal local strain of any element is 70, the strain 
amplitude. Thus a linear regime in the SGR model is ensured whenever 70 <C 1. 
The linear viscoelastic spectrum, as defined in (|l2j), can be found for the SGR 
model using (|33|): 

G*(u, t, t s ) = 1 - f e-W-^Y^Gpit - t')dt' (41) 



In principle, this quantity depends on t s , the time when the oscillatory strain 
was started. However, when the experimental timescales become large, we find 
(as shown in App. ^) that this dependence on t s is weak. In fact, within the 
SGR model, the conditions needed to make G* negligibly dependent on t s (for 
low frequencies, u<l) are that u(t — t s ) 3> 1 and ut s S> 1. The first signifies 
merely that many cycles of oscillatory strain are performed before the stress is 
measured; the second ensures that transient contributions from the initial sample 
preparation stage (the quench at t = 0) are negligible. Notably, these criteria 
do not depend on the noise temperature x, and therefore hold even in the glass 
phase (x < 1); see Fig. ^. The fact that they are sufficient even in the glass phase 
is far from obvious physically, and requires a careful discussion: we give this in 
App. [B|. Broadly speaking, these criteria are satisfied in any experiment that 
would reliably measure a conventional G*(cj) spectrum for systems with TTI. 

For the purposes of such experiments, we can therefore drop the t s argument 
and define a time-dependent spectrum G*(u, t). Our results for the long-time 
behaviour (t S> 1) of this quantity are as follows (see App. |A.3[) : 



G*(u,t) = r(x)r(2 - x^iuf' 1 for 1< x < 2 
G*(u,t) = 1 + ^- for x = l 



Ht) ' (42) 

1 



G*(u,t) = 1 - ^(iuitY- 1 for x<\ 

r (x) 



For comparison with experimental results, the simple interpolating form 

G*(ut)-i r(x)r(2 - X ) M'" 1 - 1 

may be useful; we have checked that it provides a good fit to our numerical data, 
at least over the noise temperature range to approximately 1.3 (see Fig. ^). 
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Figure 6: Demonstration of t s -independence G* (uj, t, t s ) in the glass phase. Shown 
are 1 — G'(u, t, t s ) (solid lines) and G"(ui, t, t s ) (dotted lines) against ut. The noise 
temperature is x = 0.7 and the frequency uj = 0.01. Start-time values t s obey 
Ldt s = 2 4 , 2 5 , . . . 2 9 . When u>(t — t s ) S> 1 (a sufficient number of oscillations after 
the beginning of each dataset) and wt s 3> 1 (datasets beginning further on the 
right), G* becomes independent of t s . 



By measuring G*(u>,t) we are directly probing the properties of the system 
at the time of measurement, t. In light of this, the results of (42) are easily 
understood. In the ergodic phase (x > 1), G*(u,t) will reach a t-independent 
value within a time of 0(1 /uj) after the quench, as the relevant traps will then 
have attained their equilibrium population. The relaxation time is then of 0(tq) 
(that is, 0(1) in our units) and the response G*(u,t) is a function only of uj. In 
contrast, below the glass point the characteristic relaxation time at the epoch 
of measurement is of order t, and the response is a function only of the product 
u>t. Since the losses in an oscillatory measurement arise from traps with lifetimes 
less than about 1/uj (elements in deeper traps respond elastically) , the overall 
response becomes more elastic as the system ages into traps with r > 1/uj. 

Numerical results for the viscoelastic spectrum G*(ui,t) at various measure- 
ment times t for various x are shown in Fig.|8[ These indeed show a characteristic 
"hardening" of the glassy material as it ages: the storage modulus at low fre- 
quencies evolves upwards, and the loss modulus downward (Sollich et al., 1997 



Sollich, 1998|) . Each spectrum terminates at frequencies of order uit o± 1. This 
is because one cannot measure a true oscillatory response for periods beyond the 
age of the systemPJ. Therefore, the rise at low frequencies in G' spectra like Fig. |l| 

23 And, although ( fll| ) still provides an unambiguous definition of G* (cu,t,t s ), this ceases to be 
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Figure 7: Approximate curves for G'(u;,t) and G"(u),t) generated using the sim- 
ple interpolating formula (|43| ) (thin lines) compared to the numerical data for 
these quantities (thick lines), plotted as a function of co at fixed t = 10 . Noise 
temperatures shown. Note that the predictions of the interpolating for- 

mula are practically indistinguishable from the numerical data over the frequency 
window shown. 



represents the ultimate rheological behavior^. 

It is shown in App. [B| that the insensitivity of G*(u>,t,t s ) to t s in practical 
measurements of the viscoelastic spectrum (where an oscillatory strain is main- 
tained over many cycles) arises because (even when x < 1) the most recently 
executed strain cycles dominate the stress response at time t. In essence, the re- 
sult means that, as long as oscillatory strain was started many cycles ago, there is 
no memory of when it was switched on; accordingly (by linearity) an oscillatory 
strain started in the distant past and then switched off at t s , will leave a stress 
that decays on a timescale comparable to the period of the oscillation. This is 
markedly different from non-oscillatory stresses, where long term memory implies 
that the response to a step strain, applied for a long time, persists for a similarly 
long time after it is removed (see Sec. [3| above). Thus the fact that the SGR 
glass "forgets" the t s argument of G*(uj,t,t s ), is directly linked to the oscillatory 
nature of the perturbation. As also shown in App. [B|, this forgetfulness means 
that, in the SGR model, a Fourier relationship between oscillatory and step strain 
responses is recovered; to a good approximation, one has the relation ( |i~4| ) 

poo 

G*(lj, t)=iuj e~ twt G{t', t) dt' (44) 
Jo 

Apart from the explicit dependence on the measurement timep^l t, this is the 
usual (TTI) result. But here the result is nontrivial due to the presence of ageing 

independent of t$ in this regime, so G*(u),i) is undefined. 

24 Note that this only applies for fj, = 1 in Struik's scheme, as exemplified by SGR. Whenever 
jti < 1, the region to the left of the loss peak can, in principle, be accessed eventually. 

25 Formally, t appears as the time at which an step strain was initiated, or an oscillatory 
measurement ended. Thus G*(w,i) is to within iui, the Fourier transform of the step strain 
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Figure 8: Left column: viscoelastic spectra G'(lo) (solid lines) and G"(uj) (dashed 
lines) versus frequency, u> (for x > 1) or scaled frequency u;i (for x < 1). Right 
column: frequency-dependent corrections to Hookean elasticity, 1 — G' (solid 
lines), G" (dashed lines). Data are shown for systems aged t = 10 , 10 8 , . . . 10 11 . 
At any fixed u the curves lie in order of age; data on the oldest system is marked 
by the symbols. There is good data collapse both above and below the glass 
point (with the appropriate scalings); the data for x = 1 do not collapse in either 
representation due to logarithmic terms. If plotted against u rather than ut, the 
data for x = 0.7 would resemble Fig. |]. 



31 



effects. As discussed at the end of Sec. ||, we speculate that the relation ( B4j ) 
holds not only for the SGR model, but in fact for all systems which have only 
weak long term memory. 



6.1.3 Startup of Steady Shear 

Consider now a startup experiment in which a steady shear of rate 7 1 is 
commenced at time t w . So long as we restrict attention to times short enough 
that the total strain remains small (7(2 — i w ) <C 1) the system remains in a linear 
response regime -P| 

Within the regime of linear response, any element's lifetime is independent of 
strain and obeys r = exp(E/x). As described in Sec. |5.2| above, at a time t after 
a deep quench, the distribution of lifetimes obeys P(t, t) ~ Tp(r) for r <C t and 
P(t, t) ~ tp(r) for t 3> t. Since the local stress associated with a given trap is 
of order 7T for r <C t — t w , and j(t — t w ) for r 3> t — t w , we can construct an 
estimate of the macroscopic stress; for t — t w <C t w , 



a(t) 



7 



T 2 p(r)dr + (f _ t w ) j* r/3 ( T )d T + ( t _ tw)t /«, p(T)d 



/1 Tp(r)dT + tj t °° p{r)dT 
7 [x(t - t w ) 2 - x + (x-2)(t- t^t 1 -* + x(l - x)] 



2) ft 1 



(45) 



This gives, for long times and in the linear response regime, a(t) ~ y(t — t w ) for 
x < 1 (which is purely elastic behaviour), a(t) ~ j(t — t w ) 2 ~ x for 1 < x < 2 
(which is an anomalous power law), and c(x) ~ 7 for x > 2; repeating the 
same calculation with t S> t w gives the same asymptotic scaling in each case. An 
asymptotic analysis of the constitutive equations confirms these scalings, with the 
prefactors as summarized in table |2|. Because the results depend only on t—t w , any 
explicit dependence on t w (ageing, or anomalously slow transients) must reside in 
subdominant corrections to these leading asymptotes. Accordingly, linear shear 



startup is not a good experimental test of such effects (but see Sec. 3.2.2| below). 

The power law anomaly for 1 < x < 2 can be understood by examining which 
traps make dominant contributions to a(t) = J s(r,t)dT. (Recall that s(r,t)dT 
is the stress contribution at time t from elements of lifetime r; see Sec. |4.1| .) For 
x > 2, s(r,t) is weighted strongly toward traps of lifetime 0(1); hence a(t) tends 
to a finite limit (of order 7) as t — > 00, and the viscosity of the system is finite. For 
x < 2, on the other hand, most of the weight in the s(r,t) distribution involves 
lifetimes of order t. As time passes, stress is carried by deeper and deeper traps, 
and (in the absence of flow-induced yielding) the mean stress diverges as t — > 00. 

In fact, as discussed in Sec. above, just as the Boltzmann distribution for 
the relaxation times Peq(r) = P(t, 00) ~ T p{ T ) is non-normalisable for x < 1 
(giving glassiness and ageing), so, in the absence of strain- induced yielding, is 
the ultimate distribution s(t, 00) ~ T 2 p{r) of stresses residing in traps of lifetime 

response function G(At, t) that would be measured if a step strain were applied immediately 
after the oscillatory measurement had been done. 
26 This contrasts with the ultimate steady-state behaviour which, for x < 2, is always nonlinear; 



the crossover to a nonlinear regime at late times is discussed in Sec. 3.2.2 below 
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: tor t - t w <C i w 

7 


0"(i-i w ,t w ) 

: tor t -t w ^> t w 

7 


2 < x 


x — 1 
x-2 


x — 1 
x-2 


1< x < 2 


r(x) . .n „ 
2 — x 


(x-l)r(x)^ 4 ^ 
2 - x 1 wj 


x < 1 


(t - tw) 


(l-x)(t-tw) 



Table 2: Stress response to shear strain of constant rate 7 at short and long times 
(t—t w >l,i>l,7<Cl assumed). These results apply to the regime j(t—t w ) <C 
1, where strain- induced yielding can be neglected, making the response linear. 
r(x) denotes the Gamma function. 



r, whenever x < 2. The zero shear viscosity ij is therefore infinite throughout 
this regime, as noted previously. 

6.2 Nonlinear Response 

We now turn to the nonlinear behaviour of the SGR model under imposed strain, 
starting with the step strain case. 

6.2.1 Step Strain 

The nonlinear step strain response function was defined in (0) . It is found for the 



SGR model from fl33|): 



G(t - t w , i w ; 7o ) = G (Z(t, 0)) + / Y(t')G p (Z(t, t'))dt' (46) 

Jo 



where, using (|26|): 



Z(t, t') = (t- t w ) exp ( 7 o/2x) + (i w - (47) 

On the other hand, in the linear regime we have: 

G(t - t w ,i w ,7o — > 0) = G(t-t w ,t w ) 

= G [(t - t w ) + (t w - 0)} 

+ / Y{t')G p [(t - t w ) + (t w - 0] (48) 

JO 
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Direct comparison of @ and (|4|) reveals that: 

G(t - t w , i w ; 70) = G ((t - t w ) exp (-yg/2a;) , t w ) (49) 



This result generalizes that of Sollich (1998| ) for the non-ageing case (x > 1). It 



can be understood as follows. Within the SGR model, instantaneous response 
to a step strain at t w is always elastic (that is, G(0,t w ,7o) = 1); the fraction of 
stress remaining at time t > t w is the fraction of elements which have survived 
from t w to t without yielding (see Sec. |6.1.1| above). The stress decay is therefore 
determined entirely by the distribution of relaxation times in the system just 
after the strain is applied at time t w . The effect of a finite strain is solely to 
modify the distribution of barrier heights, and hence to modify this distribution 
of relaxation times r; in fact (within the model) nonlinear strain reduces the 
yield time of every element by an identical factor of exp(7o/2x) ( tBollich, 1998 ). 



Thus the relaxation after a nonlinear step strain at i w is found from the linear 
case by rescaling the time interval t — t w using this same factor. Accordingly, the 
asymptotic results given for G(t — t w ,t w ) in table |] can be converted to those for 
the nonlinear regime by replacing the time interval t — t w by a strain-enhanced 
value (t — i w ) exp(7g/2x), wherever it appears there. 

6.2.2 Startup of Steady Shear 



In Sec. |6. 1.3 we discussed the response to start up of steady shear (with 7 <C 1) 



at time t w ; we assumed there that a linear response was maintained. Let us 
now consider the effect of strain- induced yield events, which cause nonlinearity. 
Consider first what happens for x > 2 (where the SGR model predicts Newtonian 
fluid behaviour for 7 < 1). Here the main stress contribution is from elements 
which, were they unstrained, would have lifetime t(E) = ex.p(E/x) of order unity. 
So, if the strain rate obeys 7 < 1, these elements will acquire only negligible stress 
before they yield spontaneously. Hence their lifetimes are not affected by strain, 
and the stress response remains linear at all times, including the steady state 
limit: a(t — > 00) — > 777. 

In the following, we focus on the case x < 2, where nonlinearities do appear. 
The dominant stress contributions in this noise temperature regime are from 
deep traps, i.e., elements with lifetimes of order t. Linearity applies only if 
such elements are unlikely to undergo strain-induced yielding before they yield 
spontaneously, after a time of order t. Such elements carry strains of order 7*, 
which enhances their yield rate by a factor exp[(7i) 2 /2;c]; we require that this is 
small, which holds only so long as jt <C 1. Hence the predictions of the linear 



theory of Sec. 6.1.3 can be maintained to arbitrarily long times only by taking 
the limit 7 — > before one takes the steady state limit of t — > 00. This means 
that the width of the linear response regime in steady flow is vanishingly small 
for x < 2, as previously discussed. 



As mentioned in Sec. p.l.3| , throughout the linear period the startup curve 
shows no strong ageing or transient effects, even though the stress is evolving into 
deeper traps. At finite 7, the linear period ends at t ~ 7 _1 (within logarithmic 
terms, discussed below); at later times, the main stress-bearing elements will, 
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during their lifetimes, become strongly strained. Indeed, at strain rate 7, an 
element with yield energy E will be strained to the top of its yield barrier in 
a time ti nt ~ E 1 / 2 /j ~ (logr) 1 / 2 /^. The tendency of the stress distribution 
s(r, t) (and also, for any x < 1, the lifetime distribution P(r,t)) to evolve toward 
deeper and deeper traps is thereby interrupted: the lifetime of a deep trap is 
converted from r to a much smaller value, of order (logr) 1 / 2 /^ ( |Sollich et al., 
1997i~Sollich, 1998|) . This truncation of the lifetime distribution is enough to 



ensure that these distributions are never dominated by the deep traps, and a 
steady state is recovered; accordingly, there are no ageing effects at late enough 
times either. 

Note, however, that the stress at the end of the linear regime can be higher 



than the steady state value, leading to an overshoot in the startup curve; see 3ol- 



lich (1998 ). This overshoot region, unlike the two asymptotes, shows a significant 
dependence on the system age t w , as shown in Fig. |9[ The physics of this is clear: 
the extent of the linear regime progressively gets larger as t w is increased, because 
the system has aged into deeper traps (and because the SGR model assumes that 
within each trap the relation between stress and strain is linear) . Thus the strain 
at which strong yielding sets in increases (roughly logarithmically) with t w ; the 
height of the overshoot is accordingly increased before dropping onto the same, 
t w -independent, steady-shear plateau. 








Figure 9: Stress response cr, in shear startup, vs strain 7 at noise temperature 
x = 0.3 and strain rate 7 = 0.001. Curves from bottom to top correspond to 
increasing ages t w = 10 2 , 10 3 . . . 10 9 at time of startup. 
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7 Rheological Ageing: Imposed Stress 



We now analyse the SGR model's predictions for various stress-controlled rheo- 
logical experiments. (We continue to assume the sample to have been prepared 
at time t = by the idealized "deep quench" procedure defined in Sec. |5.1| ) 
As previously remarked, the structure of the constitutive equations makes the 
analysis more difficult for imposed stress than for imposed strain. The following 
discussion is therefore largely based on our numerical results, with asymptotic 



analysis of a few limiting cases. Our numerical method is outlined in App. C.2. 



7.1 Linear Response 
7.1.1 Step Stress 

The SGR model predicts that upon the application of a step stress there will 
be an instantaneously elastic response. Elements then progressively yield and 
reset their local stresses to zero; thus we must apply progressively more strain 
to maintain the macroscopic stress at a constant value. In this way strain is 
expected to increase with time (but at a rate that could tend to zero at long 
times). Potentially therefore, individual elements can acquire large local strains 
and, just as in the shear startup case, linearity of the response need not be 
maintained at late times. As we did for shear startup, we therefore first proceed 
by assuming that the response is linear; we find the corresponding 7(i) and then 
(in Sec. [7.2| below) consider a posteriori up to what time t the linear results 
remain valid. 

In the linear regime the step stress response is described by the creep com- 



pliance J(t — t w ,t w ) which was defined for non TTI systems in Sec. 2.5. We 
computed this quantity numerically from the linearized form of the constitutive 
equation ( p3| ) for the SGR model, which for step stress may be written 

1 = J(t - t w , t w ) - f J(t' - t w , t w )Y{t')G p {t - t')dt' (50) 
Jt v 

In analysing our numerical results we first identify, as usual, regimes of short and 
long time interval between stress onset and measurement, t—t w <C t w and t—t w ^> 
t w respectively.^] In these two regimes we find the time dependences summarized 
in table |J| For the long time interval regime (t — t w 3> t w ), the results were in 
fact obtained as follows. Curves for J(t — t w , t w ) were first generated numerically; 
the observed scalings (for example, J ~ (t — t w ) x_1 for 1 < 2 < 2) were then 
taken as ansatze for analytic substitution into the constitutive equation (|50|). In 
each case this allowed us to confirm the given functional form, and to compute 
exactly the x-dependent prefactors shown. These prefactors were cross-checked 
by comparison with the numerical results; no discrepancies were found within 
available accuracy. 

To obtain results for short time intervals, we proceeded by assuming that 
the resulting compliance J(t — t w ,t w ) is the same as if we first let t w — > 00 

27 As before, we apply the "macroscopic time" conditions t — t v ^> 1 and t w S> 1- 
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J(i — t w , i w ) for t — t w t w 


J(t — t w , t w ) for t — t w » i w 


x > 2 


x - 2 . 
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(t tw) 


1< x < 2 


(t - tw)"" 1 


(« - tw)"" 1 


J- 1 €</ 11 1 i Ju I 
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x = 1 
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Table 3: Linear creep compliance in the SGR model at long and short times 
(t — t w ^> 1, i w ^> 1 assumed). T(x) denotes the Gamma function, and ^>(x) = 
r'(x)/r(x) the usual psi function. The blank entries for x < 1 are explained in 



the text; see also footnote 28. 



(the dominant traps are in Boltzmann equilibrium; see Fig. |a); this limits^ the 
analysis to x > 1. The resulting prediction of J(t — t w ,t w — > oo) was found 
analytically from G(t — t w ,t w — ► oo) and the reciprocal relations between the 
corresponding Fourier transforms (see Sec. 2.6 above); these were again checked 
numerically. 

Further insight into the results of table |3| can be gained as follows. In step 
stress, we need to keep applying larger and larger strains because elements pro- 
gressively yield and reset their local stresses to zero. To maintain constant stress, 
the rate at which stress increases due to straining, which in our units is just the 
strain rate 7, must match the rate at which stress is lost, due to local yielding 
events. The latter defines a "stress-weighted hopping rate" Y s = j 7 s(r,t)d,T. 
For x > 2, Y s remains a constant of order <7o; stress remains in traps of lifetime 
r = 0(1) and the creep response is purely viscous. For x < 2, however, Y s de- 
cays as a power law^ of (t — i w ); the stress distribution s(r,t) is dominated at 
time t by traps with lifetimes r of order t — t w , the time interval since the stress 
application. 

For 1 < x < 2, the scenario given above for the time-dependence of Y s is 
closely analogous to that given in Sec. 5^2 above for the hopping rate Y = 
J T~ l P(r,t)dT in systems with x < 1. Indeed, the evolution of Y s following a 

28 For x < 1, we find instead J = 1 + const x \{t — t w ) /t v ,] 1 ~ x at very early times; but this 
breaks down as soon as the second term becomes comparable to the leading (elastic) result. 



_>,] 



In fact, Y s ~ 7 ~ (t — t w ) v where y = x — 2 for 1 < x < 2 and y = —1 for x < 1. 
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Figure 10: Linear creep compliance J(t — t w ,t w ) against scaled time interval 
(t — t w )/t w for noise temperature x = 0.7. Curves from bottom to top correspond 
to increasing times t w = 10 2 , 10 3 . . . 10 6 of stress onset. Note the approach to 
a limiting scaling form as t w becomes very large compared with the microscopic 
time To = 1. 

step stress, at noise temperature x, is closely related^] to that of Y, following a 
quench, at noise temperature x — 1. 

The ageing behaviour of the linear creep compliance J{t — t w ,t w ) shows 
significant differences from the step strain modulus G(t — t w ,t w ) discussed in 
Sec. |6.1.1| above0. In the glass phase (x < 1), the strain response to step 
stress indeed depends on age: it is a function of (t — t w )/t w as expected (see 
Fig. 10). However, the dependence (for long time intervals) is only logarithmic; 



J(t — t w ,t w ) ~ hr((i — t w )/t w ) = ln(t — t w ) — lnt w (see table y) which means 

30 More generally one can show for the SGR model that, for an equilibrium system whose noise 
temperature is x > 1, the evolution of the stress distribution s(r, t) following application of a 
step stress at t = is, at long times, equivalent to that of the probability distribution P(r, t), in 
a system deep-quenched to a noise temperature x — 1 at t = 0. This result is connected with the 



discussion made in Sec. 3.3 above, of the variation with x of the dynamics of successive moments 
of the lifetime distribution: at noise temperature x + n, the dynamics of the nth moment is like 
that of the zeroth moment at noise temperature x. 

31 Since it refers to a shear measurement, one would not expect our result to resemble the 



empirical (stretched exponential) form measured by struik (1978) for a wide range of materials 
in tensile creep; nonetheless, it shows upward curvature on a log-log plot before approaching the 
eventual logarithmic form (with downward curvature). The same applies in nonlinear creep; see 
Fig. below. 
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Figure 11: Linear creep compliance J(t — t w ,t w ) against time interval t — t w for 
noise temperature x = 1.3 for t w = 10 3 , 10 4 . . . 10 7 (solid lines, top to bottom). 
Lower dashed line: theoretical prediction for the short time regime t — t w t w . 
Upper dashed line: asymptote for long time regime t — t w ^> t w . Note that 
short time and long time behaviours are each independent of t w (as expected for 
x > 1), but that the crossover time between them scales with i w . Inset: same 
data plotted against scaled time (t — t w )/t w ; the order of the curves is reversed. 
The crossover between short and long time behaviour takes place at a value of 
(t — t w )/i w which is roughly the same for all curves, demonstrating the scaling 
with t w . 



that in the long time interval limit (t — t w 3> t w ) the explicit waiting time depen- 
dence (lni w ) represents formally a "small" correction to the leading behaviour 
ln(t — t w ). This relatively slight ^-dependence in creep measurements is intu- 
itively reasonable: the strain response at time t to step stress is not determined 
purely by the relaxation spectrum at t w (as was the case in step strain, table [j]), 
but by the dynamics of the system over the entire interval between t w and t. This 
decreases the sensitivity to the time i w at which the perturbation was switched 



on. Similar remarks hold above the glass point (1 < x < 2, see Fig. 11): in step 



strain, we found for t — t w ^> t w a slow transient behaviour which depended to 
leading order upon t w (table |l]). For step stress, however, the corresponding i w 
dependence is demoted to lower order, and the late-time response is dominated 
by TTI terms. 

32 We restate here why we call these effects for x > 1 transient behaviour rather than ageing. 
As explained after eq. ([^J, a consistent definition of long term memory and ageing for the 
step stress response function J(t — t„,t„) requires a form of "regularization" by considering the 
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7.1.2 Oscillatory Stress 



For the SGR model it was noted in Sec. 6,1.2 that (i) in the oscillatory stress 
response G*(u,t,t s ), the t s dependence is negligible for low frequencies (lu <C 1) 
whenever uj(t — t s ) ^> 1 and ut s 3> 1; (it) these conditions are satisfied in most 
conventional rheometrical measurements of the viscoelastic spectrum, where an 
oscillatory strain is maintained for many cycles; and (Hi), perhaps surprisingly, 
these facts are true even in the glass phase, x < 1, of the SGR model. We also 
noted that, because response to oscillatory strain is dominated by memory of 
the few most recent cycles (over which the system has barely aged), G*(u,t) is 
the Fourier transform (with respect to the time interval At) of the step strain 
response function G(At, t) that would be measured if a step strain were applied 
immediately after the oscillatory measurement had been done (44). 



We have confirmed numerically that similar remarks apply to the oscillatory 
stress response function J*(u, t, i s )-0 This was defined in Sec. 2/7 as the strain 



response, measured at t, to an oscillatory stress initiated at time t s . Memory of 
the startup time t s is indeed small in J* (to, t, t s ) so long as u(t — t s ) S> 1, ujt s S> 1 
(and u<l). It appears that, just as in the case of a strain controlled experiment, 
strain response to oscillatory stress is dominated by memory to the most recent 
cycles, over which the system has barely aged. We may therefore suppress the t s 
parameter, defining a compliance spectrum at time t by J*(u,t). Furthermore, 
J*(u,t) is found numerically to be the reciprocal of G*(uj,t), 

J*(oj,t)G*(uj,t) = 1 (51) 

just as it is (without the t argument) in normal TTI systems. The numerical 



confirmation of this result is presented in Fig. |12. We emphasize that this result, 
like the previous one, has been confirmed here specifically for the SGR model; but 
it may hold more widely for systems with weak long term memory (see Sec. pf). 



7.2 Nonlinear Response 
7.2.1 Step Stress 

In Sec. |7.1,1| we argued that a step stress, o~(t) = ao@(t — i w ), of size <7o *C 1, 
induces a strain response "f(t) which increases over time, but remains linear in 
a"o for at least as long as the linearized constitutive equations predict ^(t) <C 1. 
This is because ^y(t) provides an upper bound on the local strain of each element. 
Although sufficient to ensure linearity, this is not always necessary; we require 
only that the characteristic strain of those elements which dominate the stress is 
small. For x > 2 (the Newtonian regime) the dominant elements have lifetimes 
0(1) and so the response is linear to indefinite times so long as do <C 1 (ensuring 

material in question in parallel with a spring of infinitesimal modulus g. This effectively puts 
an upper limit of J max = 1/g on the observable values of J(t — t w ,t w ). Taking the limit t w — > oo 
for x > 1 then results in a fully TTI step stress response, whatever the value of J ma x- On the 
other hand, for x < 1, the (albeit weak, logarithmic) t w -dependence of the response remains 
visible even for finite values of J < J max - 

33 Although unsurprising, this does require explicit confirmation since, for example, the tran- 
sient effects from switching on the perturbation could be different in the two cases. 
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Figure 12: Real and imaginary parts of the product G*(to,t)J*(uj,t) vs ut at 
noise temperature x = 0.7 and frequency u> = 0.01. The usual reciprocity relation 
between G* and J* is seen to hold to within about one percent. Shown are the 
results of several runs, each over a different time window. A vertical arrow marks 
the horizontal co-ordinate of the start of each data set. In each run shearing 
was commenced 20 cycles before the start of data output, to ensure that u(t — 
t s ) S> 1 (necessary for ^-independence). The oscillatory deviations, visible for the 
leftmost data set, arise because the other condition for ^-independence (tot s S> 1) 
is only just satisfied. 



■y(t) <C 1 for all times t). But, whenever x < 2, the linear analysis of Sec. [T.l.l| 
indicates the dominant elements have lifetimes of order t—t w ; so a self-consistently 
linear response is maintained only provided that j(t)(t — t w ) <C 1, just as in 
startup of steady shear (see Sec. |6.2.2 ; here we make the additional assumption 
that 7 only changes negligibly between t w and t). Using the forms for J(t — t w , t w ) 
as summarized in table ||, we then find that for 1 < x < 2 the strain response to 
step stress remains linear only for as long as t — t w <C (I/ctq) 1 ' ( x ~ x \ Beyond this 
time we expect strain-induced yielding to become important. 

To confirm the predicted linearity at short times, and to extract the long time 
non- linear behaviour, we numerically solved the non-linear constitutive equa- 
tions ( p4| , |25| ) by an iterative method (see App. G.2 ); this was done first for 
1 < x < 2 (Fig. 13). The results show a linear regime of the expected tempo- 
ral extent, followed by a crossover into a non-linear steady-state flow regime, in 
which 7(i) cx Oq t. The latter is in agreement with the flow curve (35). 

The same numerical procedure was then used for the glass phase, x < 1, for 
which the flow curve shows a finite yield stress, a y (x). As expected, the numer- 
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Figure 13: Nonlinear creep compliance J(t—t w , t w , <7o) as a function of time inter- 
val t — t w , for a step stress of size do applied at time t w = 100. The noise tempera- 
ture is x = 1.3. Solid lines, bottom to top: a = 10~ 3 , 10~ 2 - 5 , 10~ 2 , 10~ L5 , 10 -1 . 
Over the time intervals shown, the curve for uo = 10~ 3 is indistinguishable from 
the linear compliance (not shown). Dotted line: final flow behaviour predicted 
from steady state flow curve for o"o = 10~ L5 . 

ical results for step stress of very small amplitude <jq <C cr y show no crossover 
to a steady flow regime at late times. Instead, the system continues to creep 
logarithmically, according to the linear creep result (table |3|) : 

7 (t) = a J(t - t w , * w ) = °o m l _^ {x) ^ (52) 

The resulting value of j(t)(t — t w ) never becomes large; so this is self-consistent. 

Next we studied numerically the case where <jq was not small but remained 
less than the yield stress a y . For stresses not too close to the yield stress, we 
found that the creep was still logarithmic to a good approximation, but now with 
a nonlinear dependence of its amplitude on stress: j(t) rs aoA(ao)J(t — t w ,t w ). 
The prefactor A(ao) increases rapidly as do approaches the yield stress a y from 
below. Very close to the yield stress, the creep ceases to be logarithmic; ^(t) then 
grows more quickly, but with a strain rate that still decreases to zero at long 
times. On the basis of these observations, we suspect that for a given stress a"o 
the creep will be logarithmic for short times (where "short times" might mean the 
whole time window which is accessible numerically), but will gradually deviate 
from this for longer times. The deviation is expected to be noticeable sooner for 
stress values closer to yield. We attempted to verify this conjecture numerically, 
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Figure 14: Nonlinear creep compliance J(t — t w ,t w ,ao) as a function of scaled 
time interval (t — t w )/t w , for a step stress of size <7o applied at time t w . The noise 
temperature is x = 0.3. Solid curves, bottom to top: ao/cr y = 0.2 ,0.4, . . . , 1.2, 
all for t w = 100. The case do = o y is shown in bold; the dotted curve is the linear 
response result (<jq — > 0). The dot-dashed curve shows the effect of decreasing the 
waiting time to t w = 50, for ao/a y = 1.2. Comparison of the curves for the two 
different waiting times for this stress value shows that before the crossover into 
flow, the response scales with (t — t w )/t w ; once ergodicity has been restored and 
the system flows, on the other hand, scaling with t — t w is recovered. The dashed 
lines are the predictions for final flow behaviour (for the stress above yield) from 
the steady state flow curve. 



but were unable to access a large enough range of values of In ((t — t w )/t w ) to 
do so. Note that, for any cto < a y , the system ages indefinitely, and there is no 
approach to a regime of steady flow. 

Finally, as expected from the flow curve, only for stress amplitudes exceeding 
the yield stress a y (which of course depends on x) did we see an eventual crossover 
from logarithmic creep to steady flow at long times; when that happened, we 
recovered numerically the flow-curve result, 7 (t) a {a -a y ) l /^' x ^ (t-t w ). Fig.|I| 
shows examples of our numerical results that illustrate the various features of 
nonlinear creep in the glass phase mentioned above .Pj 



Note that whereas for most other shear scenarios we chose to present glass phase results for 
a noise temperature x = 0.7, we here chose x — 0.3. The yield stress is larger at this value of x, 
giving us a larger window < cto < o"y over which we see ageing and creep uninterrupted by a 
crossover into flow. 
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8 Conclusion 



In this paper we studied theoretically the role of ageing in the rheology of soft 
materials. We first provided, in Sec. [2] a general formulation of the linear and 
nonlinear rheological response functions suited to samples that show ageing, in 
which time translation invariance of material properties is lost. (Our analysis 
extends and, we hope, clarifies that of Struik (1978[) .) This was followed in Sec. || 



by a review of the concept of ageing, formally defined by the presence of long term 
memory, which can be either weak or strong. We suggested that for many rheo- 
logical applications the main interest is in systems with weak long term memory: 
these have properties that are age-dependent, but not influenced by perturba- 
tions of finite duration that occurred in the distant past. We conjectured that 
weak long term memory is sufficient to cause the age-dependent linear viscoelas- 
tic modulus to become independent of the start time t s of the oscillatory shear 
(G*(uj,t,t s ) — > G*(u>,t)) while retaining a dependence on system age t; for it to 
then obey the usual Fourier relation with the linear step strain response (likewise 
dependent on age t w ); and for it to obey a reciprocal relation G*{oj, t)J*(uj, t) = 1 
with the time-varying compliance, similarly defined. Pending a general proof of 
these conjectures, all such relationships between age-dependent rheological quan- 
tities do however require empirical verification for each experimental system, or 
theoretical model, that one studies. 

Within this conceptual framework, we then explored rheological ageing effects 
in detail for the SGR model. After reviewing the basic rheological definition 
of the model in Sec. [|, we discussed in Sec. | its ageing properties from the 
point of view of the mean jump rate Y(t) whose behaviour is radically different 
in the glass phase (noise temperature x < 1) from that in the normal phase 
(x > 1). The glass phase of the SGR model is characterized by "weak ergodicity 
breaking" , which means that the elastic elements that it describes evolve forever 
towards higher yield thresholds (deeper traps) , causing a progression toward more 
elastic and less lossy behaviour. Within the glass phase, there is a yield stress 
cr y , and for applied stresses less than this, genuine ageing effects arise. These 
phemonena were explored in depth in Sec. |6| and Sec. for the cases of imposed 
stress and imposed strain respectively. Ageing effects are distinguished from 
otherwise similar transient phenomena (arising, for example, when x > 1) by the 
criterion that a significant part of the stress relaxation, following infinitesimal 
step strain, occurs on timescales that diverge with the age of the system at the 
time of strain application. This rheological definition appears appropriate for 
most soft materials and follows closely the definition of long-term memory in 
other areas of physics (pugliandolo and Kurchan, 1995; Bouchaud et al., 1998 



Cugliandolo and Kurchan, 1993| ). 

In the glass phase of the SGR model, the nature of the ageing is relatively 
simple; for a step strain or stress applied at time t w , both the linear stress relax- 
ation function G(t — t w ,t) and the linear creep compliance J(t — t w ,t w ) become 
functions of the scaled time interval (t — t w )/t w only. This scaling is a simple 
example of the 'time waiting-time superposition' principle postulated empirically 
by [Struik (1978; ) (in the somewhat different context of glassy polymers). The 
time-dependent viscoelastic spectra G'(uj,t) and G"(uj,t) have the characteristic 
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ageing behaviour shown in Fig. [I]: a loss modulus that rises as frequency is low- 
ered, but falls with age i, in such a way that it always remains less than G'(u,t) 
(which is almost constant by comparison). For x < 1 such spectra collapse to a 
single curve (see Fig. ^) if ojt, rather than lj, is used as the independent variable. 
Note that in more complicated systems, Eq. ( p2| ) may be required instead, to 
describe ageing on various timescales that show different divergences with the 
sample age i w . Even in simple materials, there may be an additional non-ageing 
contribution to the stress relaxation which the SGR model does not have; this 
will also interfere with the scaling collapse of both G{t — t w ,t w ) and G*{uo,t). 
We found that, in its glass phase, the SGR model has weak long term memory, 
and we confirmed numerically that the conjectured relationships, Eqs. (13,14,|l5|), 
among age-dependent linear rheological quantities indeed hold in this case. 

Significant ageing was also found for nonlinear rheological responses of the 
SGR model. For example the nonlinear step-strain relaxation follows the same 
ageing scenario as the linear one, except that all relaxation rates are speeded 
up by a single strain-dependent factor (Eq. (|4^)). This form of nonlinearity is 
a characteristic simplification of the SGR model, and would break down if the 
elastic elements in the model were not perfectly Hookean between yield events. 
Another interesting case was startup of steady shear; here there is no significant 
ageing in either the initial (elastic) or the ultimate (steady flow) regime; yet, 
as shown in Fig. [9|, the intermediate region shows an overshoot that is strongly 
dependent on sample age. For an old sample, the elastic elements have higher 
yield thresholds. The linear elastic regime therefore extends further before the 
imposed strain finally causes yielding, followed by a larger drop onto the same 
steady-shear plateau. The plateau itself is age- independent: the presence of a 
finite steady flow rate, but not a finite stress, is always enough to interrupt 
the ageing process within the SGR model. Finally we found that the nonlinear 
creep compliance (Fig. |14| ), shows interesting dependence on both the stress level 
and the age of the sample; for small stresses we found logarithmic creep (for all 
x < 1), crossing over, as the yield stress is approached, to a more rapid creep 
that nonetheless appears to have zero strain rate in the long time limit. Non- 
linear creep gives challenging computational problems in the SGR model, which 
is otherwise simple enough, as we have shown, that almost all its properties can 
be calculated either by direct asymptotic analysis or using (relatively) standard 
numerics. Remaining drawbacks include (from a phenomenological viewpoint) 
the lack of tensorial elasticity in the model and (from a fundamental one) uncer- 
tainty as to the proper physical interpretation, if one indeed exists, of the noise 
temperature x flSollich et al., 1997; Sollich, 199^ ). 

Though obviously oversimplified, the SGR model as explored in this paper 
may provide a valuable paradigm for the experimental and theoretical study 
of rheological ageing phenomena in soft solids. More generally, the conceptual 
framework we have presented, which closely follows that developed to study age- 
ing in non-flowing systems such as spin-glasses, should facilitate a quantitative 
analysis of rheological ageing phenomena across a wide range of soft materials. 
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A Calculation of Linear Response Properties 



A.l Initial Condition 

In discussing the SGR model's non-equilibrium behaviour (Sees. || and ^) we 
considered for definiteness a system prepared by a quench from an infinite noise 



temperature (see Sec. 5.1), i.e., with an initial distribution Pq(E) = p{E) of yield 
energies or trap depths. For our predictions to be easily compared to experi- 
mental data, however, they must be largely independent of the details of sample 
preparation. To test for such independence, we consider the extent to which our 
results would change if the pre-quench temperature, which we denote by xq, were 
finite. This corresponds to an initial trap depth distribution 

P (E) oc exp(E/x )p(E) (53) 

In this appendix, we restrict ourselves to the linear response regime, where the 
effects of finite xo (if any) are expected to be most pronounced; nonlinearity 
tends to eliminate memory effects. The same is true for high temperatures, and 
correspondingly we will find that the influence of xq on our results is confined 
mainly to final (post-quench) temperatures x within the glass phase (x < 1). 

A. 2 Yield Rate 

The yield or hopping rate is the basic quantity from which other linear response 



properties can be derived; see eqs. (£38,41,[5C]). It can be calculated from the 
second constitutive equation (f25|) 

1 = G (t) + ! Y(t')GJt - t')d£ (54) 
Jo 

where we have replaced Z(t, t') by t — t', as is appropriate in the linear response 
regime. The function Go(t) is defined in (p7|); for the initial condition (|53| ) it is 
related to G p via 

G (t) = G p (t,y), y = x{l-l/x ) 
where we have now included explicitly the noise temperature argument (y) in the 



argument list of G p . Substituting this into (54), and taking Laplace transforms 



with A as our reciprocal time variable, we get: 

1 



and hence 



in which (taking Laplace transforms of (|27|)) 



x G p (X,y)+Y(X)G p (X,x) (55) 



OO — X — 1 (-OO 



G p (X,x)=xl i JT - rT dr = xJ i TTJ -dr (57) 
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In its present form (56) cannot be inverted analytically. We will focus on the 
long time regime, however, where progress can be made by using an alternative 
expression for G p . From fl57j), G p (X,x) has poles at A = — r _1 . Because of the 
integration over all r = 1 . . . oo, these poles combine into a branch cut singularity 
on the (negative) real axis between A = — 1 and A = 0. We will now derive 
an expression for G p that is valid near this branch cut. This expression does 
introduce spurious singularities on the negative real axis for A < 1. But after 
inversion of the Laplace transform these only give contributions to G p (t) decaying 
at least as fast as exp(— t); they can therefore be ignored in the long-time limit. 
We first write fl57j) as 



l - r°° t~ x r 1 t 

- G p (X, x) = -^dr- — — dr (58) 
x Jo 1 + At Jo 1 + At 



After the rescaling Ar — * r, the first term becomes a representation of the Beta 
function.^] In the second term, because now r < 1, we can expand the de- 
nominator into a series that is convergent for |A| < 1. This gives the desired 
expression 

oo 

G P (X, x) = a{x)X x ~ 1 + K(x)X n (59) 



n=0 



in which 



a{x)=xF(x)T(l-x), b n (x) = ^—t (60) 

n + 1 — x 



This is valid for |A| < 1 and therefore in particular near the branch cut A = 
— 1 . . . 0; in the representation (j59|), this branch cut is apparent in the fractional 
power of A in the first term. The above derivation applies a priori only for x < 1, 
because otherwise the integrals in (|58| ) diverge at the lower end. However, using 
the relation 

-L-G p (X,x + l) = ---G p (X,x) 

X + 1 XX 



which follows directly from (|57|), it can easily be shown that fl59[ ) holds for all 
x. (For integer x, there are separate singularities in the first and second term 
of (|59D, but these just cancel each other.) 



We can now substitute ( |5^ , |60|) into (|5q) and expand the denominator to find 
a readily invertible expression for Y(X). Clearly the manner in which we perform 
the expansion depends on whether x > 1 or x < 1. Abbreviating a(x) = a, 
a(y) = a', and b n (x) = b n , we have for x > 1: 

Y{x) = x[b 6g" b 

which, upon inversion of the Laplace transform, gives: 



I _ "a- 1 - ^xy + o { a 2 ^- 1 ), \v+*-\ x, . . .) 

bn \ / 



v ; b r(2 -x)bl r(i - y) b G V ' ; v > 



35 The rescaling can be carried out only when A is real and positive. But by analytic continu- 
ation, the result (^) also holds for complex A outside the branch cut of A 2 ^ 1 , i.e., everywhere 
except on the negative real axis. 
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the first term of which is the asymptotic expression for Y(t) above the glass 
points, as in (j3(f). For x < 1 on the other hand, we have 



y(A) 



A 



l-x 



-l-x 



+ 



o(a 3 M A^+ 2 M,...) 



(62) 



which can be inverted to give 



Y(t) 



1 t 



x—l 



b t 2 ^ 



r(x) a r(l + 2(x-l)) 



aT(x - y) 



t x-i-y +0 ^(x-l) jt 2(x-l)- 
(63) 

the first term again being in agreement with (|36|). Finally, to obtain 5^(i) the 
glass point x = 1 we rewrite (|62| ) as: 



y(A) 



i 



z n W + O A, Afr+W 1 -*) . . .) 



,n=l 



in which z(A) = — 60A 1 x /a and p is the largest integer which is less than 1/(1— x) 
Inversion of the Laplace transform gives 



-1 p 



z n {t) 



+ o t~y,t 



y Jp+l)(x-l) 



bo ^r(l + n(x-l)) 

in which z(t) = —bot x /a. The Gamma function can now be expanded around 
A = 1 — x = 0; the sum over p can be performed explicitly for each term in this 
expansion. Retaining only the dominant terms for small A, and also taking the 
limit A — > of the quantities z(t), a and bo, one finds eventually 

1 rm / 1 \ 

lim Y(t) = T7 - + -i/r + O - T — 
x-i ]n(t) ln 2 (t) \]n 3 (t) J 

as stated in (j36|). 

Consider now the effect of the pre-quench temperature xo on the above re- 
sults for the asymptotic behaviour of the hopping rate Y(t). We note first that 
all the leading terms are independent of y and hence of xq. For x > 1, the 
largest y-dependent subleading term (t~ y ) in (|6l| ) becomes more important for 
smaller pre-quench temperatures xq. However, provided we restrict ourselves to 
the regime xq > x (i.e., to a non-equilibrium situation in which a quench is 
actually performed; x = xq corresponds to equilibrium conditions), we see that 
y > x — 1 and that, even to subleading order, Y{t) is independent of xq. (We 
note furthermore that in the case of the deep quench defined in Sec ^1], y = x 
and the term t~ y is very small.) For x < 1, in (p3]), the relative importance of 
the largest y-dependent term (i^ -1- ^) again depends upon the relative values 
of the pre- and post-quench temperatures. For a high enough pre-quench tem- 
perature (specifically, provided y > 1 — x, i.e., provided xq > x/(2x — 1)) the 
leading and subleading terms of Y(t) are independent of xq. For any post-quench 
temperature x < 1/2, the subleading term necessarily depends upon xq since the 
condition defined above for independence cannot be satisfied. (This is physically 
intuitively reasonable, since in general we expect a system at a lower temperature 
to remember its initial condition more strongly.) 
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A. 3 Step Strain and Oscillatory Strain Response 

Once the yield rate Y(t) is know, the linear stress response G(t — t w , t w ) to a step 
strain can be calculated from (|38|). To get its asymptotic behaviour for t — t w S> 1, 
t w S> 1, the two regimes in which the time interval t — £ w is much less and much 
greater than the age at the time of stress application t w have to be considered 
separately. In the first regime ), one can Taylor expand the hopping 

rate Y around its value at time t. In the second regime, we rewrite (|38|) as 

G(t - t w , t w ) = G (t) + / Y(t')G p (t - t')dt' (64) 

Jo 

The first term on the right-hand side can then be shown to be subdominant 
(at least for xq — > oo; see Sec. A.4| below), and the second can be treated by 



expanding G p (t — t') around t' = 0. To leading order, one then finds the results 
in table |l[ The asymptotic behaviour of the stress response to oscillatory strain, 
G*(u, t, t s ), is obtained in a similar manner from (|4l|). 



A. 4 Rheological Irrelevance of Initial Condition 



In App. [A.2| we discussed the influence of the initial state of the sample, as pa- 
rameterized by the "pre-quench" temperature xq, on the yield rate Y(t). Now 
we consider the effects of xq on the various (linear) rheological observables, con- 
centrating on the regime x < 1 where such effects are expected to be most pro- 
nounced. We begin with the response to a step strain, G(t — t w , t w ). In the short 
time regime t — t w <C t w , it follows directly from (^) that xq affects only sub- 
dominant terms (through its effect on Y(t)). In the long time regime t — t w ^>t 



w - 



we see similarly from (64) that any effect on the leading behaviour can only be 
through the first term on the right-hand side, Go(t) = G p (t, y) ~ t~ y . Comparing 
this with the second term, which from table [l| is ~ (t w /t) x (note that t ~ t — t w in 
the long time regime), and using y = x(l — 1/xq), one finds that the effect of xq is 
negligible up to t « t^°. For larger t, G(t — t w ,t w ) rs Go(t) w Go(t — t w ) and the 
response is TTI to leading order. An intuitive explanation for this behaviour can 
be found by analysing the evolution of the relaxation time distribution P(r,t w ) 



with t w ( Fielding, 2000 ). It can be shown that the initial condition P(t,0) is 
remembered in the long time tail of this distribution, t > tj. For times t>^°, 
these long relaxation times dominate the behaviour of G(t — t w , t w ) and cause the 
observed ^-dependence. 

For the step stress response J(t — t w , t w ), we found in Sec. [7.1. 1| that memory 
effects are rather weaker than for the step strain response. This is because J is 
sensitive to the average behaviour of the relaxation time distribution P(t, t') over 
the time interval t' = t w . . .t, while G depends on P(t, t w ) only. Correspondingly, 
we also find that J(t — i w , t w ) is affected only weakly by the initial preparation of 
the system and hence by xq. All effects are in subdominant terms; for the long 
time behaviour in the glass phase, for example, one finds that the asymptotic 
behaviour J(t — t w ,t w ) ~ ln((£ — t w )/t w ) is only changed by an xo-dependent 
constant offset ( |Fielding, 200C| ). 

Finally, consider the oscillatory response functions G*(u,t,t s ) and J(u,t,t s ). 
Any linear oscillatory perturbation effectively probes only those traps which have 
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a relaxation time r < 1/u. Provided such traps have attained an xcr independent 
distribution by the time the perturbation is switched on at t s , G*(u>,t,t s ) and 
J*(u>,t,t s ) will be insensitive to xq. It can be shown flFielding, 2000[ ) that the 
requirement for this is r <C t x ° for all t < 1/tu and hence ojt x ° S> 1. We argue 
in App. |B|, however, that in order to get a sensible measurement of G* (and J*) 
which is independent of start time t s , we must ensure ut s 3> 1. This condition 
then automatically guarantees that the results are independent of xq. 

In summary, the only significant effects of the initial sample preparation ap- 
pear in the step strain response at long times (t t^). In the other linear 
response properties that we studied, the initial condition only affects subdomi- 
nant terms. We reiterate our earlier statement that for nonlinear response, the 
initial sample condition should be even less important, because nonlinearities 
tend to wipe out memory effects. 



B Irrelevance of Switch-on Time in the Glass Phase 



It was stated in Sec. |6.1.2 that G*(u),t,t s ) does not depend on t s so long as 



u)(t — t s ) ^> 1 and cvt s 3> 1. These criteria do not depend on the noise temperature 
x, and therefore hold even in the glass phase, x < 1, where ageing occurs. 

This behaviour can be understood as follows. Consider a material which has 
not been strained since preparation except during a time window of duration t* 
before the present time t. First write the linearized constitutive equation as: 



, lt)= -^m^i^ ( 65) 



dG(t - 1', t' 

7^ ) 

it-t 

where, for the SGR model 

dG(t-t',t') 



dt' 



6(t-t') + Y(t')GJt-t') (66) 



with0 G p {t-t') ~ {t-t')~ x . 

Now consider the case of a step strain imposed at t—t*, so that ^(t) is constant 
in (]6q). Because dG/dt' contains a contribution of order (t— t f )~ x , the integral has 
significant contributions from t' near t — t* whenever x < 1; in fact in the absence 
of the factor Y(t') the integral would not even converge to a finite limit as t* 
becomes large. This is a signature of long-term memory: Even the strain history 
in the distant past has an effect on the stress at time t. On the other hand, for an 
oscillatory strain (likewise switched on at t — t*) one has ( |65| ) with j(t) = "foe lu3t , 
and even without the factor Y(t') the integral would now converge to a finite 
limit so long as uot* is large. The convergence of the oscillatory integral follows 
from the mathematical result known as Jordan's lemma ( popson, 1962 ) which, 



crudely speaking, states that inserting the oscillatory factor e lult has a similar 
effect to converting the integrand, dG/dt', to u>d 2 G/dt' 2 . (Physically, this extra 
time derivative arises since the stress at t due to any previously executed strain 
cycle must involve the change in dG/dt' over the cycle: if this change is small, 

36 This result follows by differentiation of respecting the fact that G(t — t„,t„) vanishes 
for negative t — t w (that is, it contains a factor Q(t — t w ) which is conventionally suppressed). 
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the response to positive and negative strains will cancel.) Accounting for this 
extra time derivative, it is simple to check that the most recently executed strain 
cycles indeed dominate the response at time t, in contrast to the non-oscillatory 
case where the entire strain history contributes. 

This observation allows us to simplify ( |65| ) further by setting G(At, t') — » 
G(At, t) where At = t — t' (we assume ut 3> 1, so that the variation in the stress 
response function over a fixed number of recent cycles is negligible). Likewise, 
the limit of integration can safely be set to A = oo. Thus we have 

G>, t) = | o °° e-- At dAt (67) 

This can be integrated by parts to give ( flij) as required. 



C Numerical Methods 

C.l Yield Rate in the Linear Regime 

To obtain numerical results for the linear response properties of the SGR model, 
the yield rate Y(t) has to be calculated first. A convenient starting point for this 
can be obtained by differentiating (| 



Y(t) = -G' (t) - J Y(t')G' p (t - t') dt' (68) 

This is a Volterra integral equation of the second kind, which can in principle be 
solved by standard numerical algorithms ( Press et al., 199^ ). Such algorithms are 



based on discretizing the time domain into a grid to = 0, t% . . . t n ; the values Y n = 
Y(t n ) are then calculated successively, starting from the known value of Yq. The 
subtlety in our case is the choice of the grid: Because for times t>lwe expect 
the hopping rate to be a power law, we expect relative discretization errors given 
a time-step At to scale roughly as At/t n . Once we have chosen an acceptable 
(constant) value for the discretization error we are therefore at liberty to increase 
the time-step At linearly with the time t n , which corresponds to using a geometric 
time grid. This allows us to generate data over many decades without too much 
computational effort. To improve accuracy, we also used a spline interpolation 
between the known points (to,Yo), (ii, Yi) . . . (i n _i, Y n _i) when determining the 
next value Y n . 



C.2 Strain Response to Finite Step Stress 

The numerical scheme used to solve (^4|) and (pq ) in the case of an imposed 
step stress a(t) = o"oG(t — t w ) is rather more complicated. This is because both 
the strain and the hopping rate, which are coupled through nonlinear integral 
equations, have to be calculated as functions of time. Again, we discretize time 
into a grid to, t\ . . . t n , where to = ^w; an d proceed along the grid calculating the 
strain 7 n and the hopping rate Y n for successive values of the index n.p*| The first 

37 Note that the integral form of the constitutive equations renders the strain and the hopping 
rate at any time step t n dependent upon the values of these quantities at all previous times 
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data point 70 = "f(t^) an d Yq = ^(*w) on this grid is then obtained directly by 
treating the discontinuity at t w "by hand" . At any subsequent time-step the two 



non- linear constitutive equations (24) and fl2q ) are solved simultaneously. The 
first is essentially of the form: 

= f(ln,Y n ,{j n ,},{Y n ,},t w ,a ) for < n' < n (69) 

while the second can be differentiated and rearranged to give 

Y Tl = g(j n ,{ln'},{ Y n'},t w ,cr ) for < ri < n (70) 



Because (|69D cannot be solved explicitly for 7„, we use an iterative process. At 
each time-step we start by placing sensible upper and lower bounds on 7„, derived 
from physical expectations about the time dependence of the strain 7(i). Each 



bound in turn is substituted into (70) (to find the corresponding value of Y n ) 
and (with its Y n ) into the function / of the right hand side of d69|). The secant 
method (press et al., 1992] ) is then used to update one of the bounds, and the 
new bound used to calculate a new Y n and /. This process is repeated until we 
obtain a sufficiently small value of / (|/| < 10~ 8 ). The current values of 7„ and 
Y n are then accepted and we proceed to the next time-step. 

We initially chose a geometric grid of time values to,ti . . ., but this led to 
numerical instabilities. We therefore switched to an adaptive procedure which 
chooses time-steps such that the strain always increases by approximately the 
same amount in a given time-step. 

Finally, note that at each iteration loop of each time-step we in principle need 
to evaluate double integrals of the form I = Jq h(Z(t,t'))dt' in which 

Z(t,t') = [ t dt"exp{[ 1 (t")- 1 (t')] 2 /2x} (71) 
Jt' 

Because this is very costly computationally, we first calculate at each loop Z(t, t') 
on a grid of t' values ranging from to t and set up an interpolation over the 
calculated points. We are then left with single integrals of the same form as I, 
and look up the value of Z whenever the integrand is called. 
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